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Abstract 

One of the most cited papers in Applied Mechanics is the work of Eshelby from 1957 
who showed that a homogeneous isotropic ellipsoidal inhomogeneity embedded in a ho¬ 
mogeneous isotropic host would feel uniform strains and stresses when uniform strains or 
stresses are applied in the far-held. Of specihc importance is the uniformity of Eshelby’s 
tensor S. Following this paper a vast literature has been generated using and developing 
Eshelby’s result and ideas, leading to some beautiful mathematics and extremely useful re¬ 
sults in a wide range of application areas. In 1961 Eshelby conjectured that for anisotropic 
materials only ellipsoidal inhomogeneities would lead to such uniform interior helds. Al¬ 
though much progress has been made since then, the quest to prove this conjecture is still 
not complete; numerous important problems remain open. Following a different approach 
to that considered by Eshelby, a closely related tensor P = SD° arises, where D° is the 
host medium compliance tensor. The tensor P is associated with Hill and is of course also 
uniform when ellipsoidal inhomogeneities are embedded in a homogeneous host phase. Two 
of the most fundamental and useful areas of applications of these tensors are in Newtonian 
potential problems such as heat conduction, electrostatics, etc. and in the vector problems 
of elastostatics. Knowledge of the Hill and Eshelby tensors permit a number of interesting 
aspects to be studied associated with inhomogeneity problems and more generally for inho¬ 
mogeneous media. Micromechanical methods established mainly over the last half-century 
have enabled bounds on and predictions of the effective properties of composite media. 
In many cases such predictions can be explicitly written down in terms of the Hill, or 
equivalently the Eshelby tensor and can be shown to provide excellent predictions in many 
cases. 

Of specific interest is that a number of important limits of the ellipsoidal inhomogeneity 
can be taken in order to be employed in predictions of the effective properties of e.g. layered 
media, fibre reinforced composites, voids and cracks to name but a few. In the main, results 
for the Hill and Eshelby tensors associated with these problems are distributed over a wide 
range of articles and books, using different notation and terminology and so it is often 
difficult to extract the necessary information for the tensor that one requires. The case 
of an anisotropic host phase is also frequently non-trivial due to the requirement of the 
associated Green’s tensor. Here this classical problem is revisited and a large number of 
results for problems that are felt to be of great utility in a wide range of disciplines are 
derived or recalled. A scaling argument leads to the derivation of the Eshelby tensor for 
potential problems where the host phase is at most orthotropic, without the requirement of 
using the anisotropic Green’s function. Goncentration tensors are derived for a wide variety 
of problems that can be used directly in the various micromechanical schemes. Both tensor 
and matrix formulations are considered and contrasted. 


1 Introduction 

The canonical isolated inhomogeneity problem has been of fundamental importance in a number 
of materials modelling problems now for well over a century. This problem is the following: a 
single inhomogeneity, i.e. a particle of general shape, with different material properties to that 
of the surrounding material is embedded inside an unbounded (in all directions, i.e. free-space) 
homogeneous host medium. Given some prescribed conditions in the far-held, what form do 
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the fields take within the inhomogeneity? As well as being interesting in its own right, this 
problem is of utmost importance in homogenization, micromechanics and multiscale modelling. 

The hrst to consider this kind of inhomogeneity problem was Poisson in 1826 |92j who 
studied the perturbed field due to an isolated ellipsoid in the context of the Newtonian potential 
problem. He showed that given a uniform electric polarization (or magnetization), the induced 
electric (or magnetic) field inside the ellipsoid is also uniform. In 1873 Maxwell [7l] derived 
explicit expressions for this held. Early work in linear elasticity saw a number of studies 
determine the held inside and around inhomogeneities, including the important case of a cavity 
(since this was correctly recognized as a defect or haw). Examples of these works were those 
associated with the case of spheres [105] , [33] , spheroids [22] and ellipsoids [101111021195] but all 
considered specihc loadings, usually of the homogeneous type in the far held, meaning uniform 
tractions or displacements that are linear in the independent Cartesian variable say x. 

The inhomogeneity problem is now usually associated with the name of Eshelby because in 
1957 he showed that for general homogeneous conditions imposed in the far held, the strain set 
up inside an isotropic homogeneous ellipsoid is uniform [23|. In 1961 Eshelby [25| conjectured 
that ...amongst closed surfaees, the ellipsoid alone has this convenient property....”. Is this 
true? In the sense of what it is thought that Eshelby meant when he made this conjecture (the 
so-called weak Eshelby conjecture, where the interior held must be uniform for any uniform 
far-held loading), this statement certainly is true although this was only proved in 2008, 
simultaneously by Kang and Milton [3B] and Liu (BS] in the case of isotropic media. There is 
a slightly different version (the so-called strong Eshelby conjecture), where the interior held 
must be uniform only for a speeific, single uniform far-held loading. This strong conjecture has 
still not been proven in the context of three dimensional isotropic linear elasticity, although 
signihcant progress has been made in the last decade, see [l5| for a review. Eurthermore the 
results obtained in [1] go beyond the weak Eshelby conjecture but still do not fully prove the 
strong conjecture. Interestingly the associated (weak) conjecture for the Newtonian potential 
problem was proved some time before Eshelby’s 1957 elastostatics paper, by Dive in 1931 m 
and Nikliborc [85] in 1932, see also the discussion in [36], [63|, [l5]. In deriving these results, 
Dive and Nikliborc proved the converse of Newton’s theorem that if V is an ellipsoid of uniform 
density, the gravitational force in V is zero |48] . The strong conjecture in the context of the 
potential problem is true in two dimensions [22] but is not true in dimensions greater than 
two. A non-ellipsoidal counterexample associated with a specific far-held loading (equivalently 
a specihc eigenstress) was found by Liu [63] ■ 

It is important to note that the proofs of Eshelby’s conjectures in elastostatics referred 
to above correspond to simply connected, isotropic inhomogeneities with Lipschitz boundaries. 
Eshelby’s work was followed up with work by numerous researchers who considered the general 
anisotropic case m, ra, [Ml, m, isni, m, la, msi- In 1974 Cherepanov [T5] 

proved that multiple inhomogeneities of non-ellipsoidal shape can interact in order to render 
the interior helds uniform; see also Kang and Milton [36] and Liu [63] who coined the term 
E-inclusions for such interacting inhomogeneities. Liu and co-workers have also considered 
the periodic Eshelby problem in two dimensions [62] . [61] . Kang and Milton [36] used their 
approach to prove Eshelby’s weak conjecture in the context of the fully anisotropic potential 
problem. Most notably, it is stressed again that the weak Eshelby conjecture for elasticity has 
not yet been proved in the context of anisotropic elasticity. 

Interest in deriving the Eshelby tensor for non-ellipsoidal inhomogeneities has always been 
present in order to show that the conjecture holds for specific classes of inhomogeneities. 
Particular attention has been paid to polygonal and polyhedral inhomogeneities and the asso¬ 
ciated properties of Eshelby’s tensor [83], [96], [81], [86], [68], [67], [63], [37]. The supersphere 
case has been considered recently by m building on the work by [87], [88], [89]. A general 
method was developed by Ru [98] in order to obtain an analytical solution associated with 
a two dimensional inhomogeneity of arbitrary cross section and explicit forms of the stress 
inside hypotrochoidal and rectangular inhomogeneities were derived. Some analytical expres¬ 
sions have recently been derived for two-dimensional problems in the Newtonian potential and 
plane elastostatics problems where inhomogeneities are either polygonal or their shape can 
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be described by finite Laurent expansions |125] . [126] , Furthermore useful properties of the 
Eshelby tensor have been deduced, including the relationship of the averaged Eshelby tensor 
for non-ellipsoidal inhomogeneities to their ellipsoidal counterparts 

More recently the inhomogeneity problem has been studied in the nonlinear elasticity con¬ 
text where in two dimensions results associated with Eshelby’s conjecture have been proved 
in two dimensions for so-called harmonic materials M, m, [50]. Although nonlinear prob¬ 
lems are generally more difficult that linear elastostatics, the nonlinearity frees up a number 
of issues that are more constrained in linear problems. The study of nonlinear problems with 
dilatational eigenstrain was recently carried out in Giordano m considered the non- 

linearly elastic inhomogeneity problem but where the constitutive behaviour is described via 
expansions in strain (Landau elasticity). 

Here attention is restricted to linear problems for ellipsoidal inhomogeneities and associated 
limits. A general approach to deriving the Hill tensor and proving many of its properties is to 
use the integral equation form of the governing equations In fact Eshelby approached the 

problem in quite a different manner, using the concept of eigenstrain [23]. Hill [37] considered 
the so-called polarization (hence P) of an ellipsoid. The review articles of Walpole [109] and 
Willis who developed the integral form of the P-tensor have been very influential and 

the text of Mura [50] describes the associated Green’s tensor and form of Eshelby tensors for 
elastostatics in detail. The consideration of isolated inhomogeneity problems allows the deriva¬ 
tion of so-called concentration tensors for dilute micromechanical schemes, where interactions 
between inhomogeneities are not important |119j . In the field of micromechanics a number 
of very ingenious approximations have been made that lead to rather excellent predictions of 
effective properties in the case where interactions amongst inhomogeneities are important (see 
e.g. [II3], [69] for broad overviews). Finally it is noted that variational bounds can be 

conveniently written down in terms of the Hill or Eshelby tensors [35], [m], [nz], [93], [TO], 

m- 

There is no real preference for the direct integral equation approach leading to the Hill 
tensor, over the Eshelby eigenstrain approach. It is chiefly down to individual preference 
although it is important to note that Hill’s tensor possesses the major symmetries whereas 
Eshelby’s does not in general. Some find the notion of eigenstrain rather artificial, although 
in many cases it is a very useful concept as a means for solving harder problems such as the 
case of multiple inhomogeneities m, CM]. The simple relation 

S = PC° (1.1) 

between the Hill (P) and Eshelby (S) tensor, where is the host modulus tensor, means that 
deriving one immediately yields the other. 

The Hill and Eshelby tensors are of great utility in a number of micromechanical methods 
and what is quite astonishing is that they can be evaluated analytically in a large number of 
very important cases. However, results are distributed over a large number of articles, reviews 
and textbooks, and furthermore often in articles that span a wide range of scientific fields due 
to the wide ranging applicability of the theory. References dealing with derivations of specific 
results are those of dz], m\, [S2], m, 0, m and [55]. The field is still very much alive, 
pushed forward by both unresolved theoretical issues as well as applications involving not only 
inhomogeneities but also cracks and dislocations |82j . [123j and by the desire to fully resolve 
the open issues described above. Recent work has focused in more detail on inhomogeneities of 
general shape and how these can feed into models of inhomogeneous media with distributions 
of non-canonical inhomogeneities 0, 0, m, CM], CIS], [126]. Such studies are important 
to understand how local stress fields develop in the medium under loading. This is highly 
dependent upon the inhomogeneity shape. 

Here the objective is to gather together important results associated with the Hill and 
Eshelby tensors for ellipsoidal inhomogeneities in consistent notation, derive a number of im¬ 
portant limiting cases such as those associated with cracks and cavities, derive compact results 
associated with the anisotropic potential problem and finally derive and state associated con¬ 
centration tensors. This should prove useful to many who frequently require the form of the 
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P- or S-tensors in practice but who struggle to find the appropriate reference. 

An important point to note is that using the so-called invariant notation, potential and 
linear elastostatics problems can be considered simultaneously, only that the latter is a higher 
order tensor analogue of the former. Here however the applications are made distinct to stress 
the different results and mechanisms for deriving these expressions. In particular the results 
from potential theory feed into those from linear elastostatics. As a result index notation shall 
be used almost entirely throughout. 

In much of the literature on micromechanics the terms inclusion and inhomogeneity are 
used interchangeably. However in some cases they are used to make an important distinction. 
An inhomogeneity is defined as a particle of general shape having different material properties 
to those of the surrounding medium in which it is embedded. On the other hand the termi¬ 
nology inclusion is used to represent a general shaped region within some medium that has 
the same properties as the surrounding medium but where this finite inclusion region has been 
subject to some eigenstrain (e.g. thermal strain). This differentiation is used e.g. in Mura [80] 
and Qu and Cherkaoui [M] . 

In ^ the integral equation formulation of the inhomogeneity problem is stated, yielding 
integral equations for the potential gradient and strain inside an inhomogeneity. In ^ it is 
illustrated that such fields are uniform when the inhomogeneity is ellipsoidal and the general 
expressions for the associated Hill tensors are stated. The notion of concentration tensors is 
also discussed. In §21 and [S] specific results are then stated and derived for the cases of the 
Newtonian potential problem and elastostatics respectively. A closing discussion is given in 
21 describing how the results are used in micromechanical methods together with a summary 
of current areas of associated research. Numerous important details and results are stated in 
Appendices in order for this review to be comprehensive but also to aid the flow of the reader. 

As many pertinent references are given as possible; the focus is specifically on the for¬ 
mulation of the Eshelby, Hill and concentration tensors rather than articles associated with 
micromechanical methods, of which there are thousands. For the latter the interested reader 
is referred to the many textbooks that have been written over the last decade, see e.g. 

mi, m. m- 


2 Integral equation formulation 

Index notation shall be used for tensors throughout, working in Cartesian coordinates and using 
repeated subscripts to imply summation. The term unbounded will be used when referring to 
free-space, i.e. unbounded in all directions. Although a general invariant formulation can be 
employed to deal with problems in the potential and linear elastostatics context simultaneously 
m , this approach can obfuscate details that are important when it comes to deriving specific 
Hill and Eshelby tensors for given anisotropies and inhomogeneity shapes. 

Notation is as defined in Fig. [1] for both the potential problem and linear elastostatics. A 
single isolated inhomogeneity Vi, for the time being of general shape and with surface dVi is 
embedded (perfectly) inside an unbounded homogeneous medium V and we denote the medium 
exterior to Vi as H\Vi = Vq. Both materials are considered generally anisotropic so that their 
material modulus tensors are 


Qj(x) = + Cfji^ - xHx)) (2.1) 

in the context of the potential problem and 

Cijkei^) = - X^x)) (2.2) 


in the context of elastostatics. Here the so-called characteristic function associated with a 
domain Vi, has been employed, being defined as 


X'(x) = 


4 


X G El, 
X ^ El. 


(2.3) 




Finally it is noted that the inhomogeneity and and host {Cfj and modulus 

tensors are uniform tensors, meaning that each component of the tensor is constant but these 
constants can be different. 



Figure 1: An inhomogeneity Vi of general shape and with boundary dVi is embedded perfectly 
inside the host medium Vq. The classical inhomogeneity problem is to determine the fields 
that arise within the inhomogeneity and host medium given some far-held condition. 


2.1 The potential problem 

Since it is often useful to consider a specihc physical problem, certainly in terms of language 
and terminology, the potential problem is described in the context of steady state thermal 
conductivity. The equation governing the steady state temperature distribution r(x) in the 
medium described above and depicted in Fig. [T] is 

d f dT\ 

where we note that no heat sources are present. The free-space Green’s function associated 
with the host phase satishes 

( 2 , 5 ) 

as well as the far-held condition limx^oo G{x.) = 0. Assuming continuity of temperature 
and normal hux across dVi, the resulting temperature distribution may be straightforwardly 
derived in integral equation form as 

f f)T BC 

T(y) = T-(y) - (Cl, - Ct,) a^(x)^(x - y) dx (2.6) 

which holds for all y. Here T*{y) is the solution to the equivalent problem satisfying (|2.4|) 
with no inhomogeneity present (or equivalently with = G? ). Upon taking derivatives of 
(12.61) with respect to yi and noting the property BGjdxi = —BG/dyi it is found that for all 
yeG, 

ei(y) = e- (y) + [Clj - ^kj) Qy,Qy, efc(x)G(x - y) dx (2.7) 

where the ith component of the temperature gradient has been dehned as e, = BT/dxi. 

2.2 Elastostatics 

The origins of the P-tensor reside in the context of elastostatics rather than in potential 
problems even though the theory is of course analogous. The P-tensor originated with Hill 
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m who also introduced the compact notation (now commonly referred to as Hill notation) 
for transversely isotropic (TI) fourth order tensors, which we summarize in Appendix IC.4.31 
Walpole [108] , Willis |115[ 111611117| and Laws [52j amongst others followed this with influential 
work associated with inhomogeneities of specific shapes, paying particular attention in many 
cases to the scenarios of discs, fibres and cracks. A number of P-tensors are also stated in 
the excellent concise review of micromechanics by Markov [69] although unfortunately, some 
typographical errors are present there and we correct those here. 

The solution to the isolated inhomogeneity problem in elastostatics proceeds analogously 
to the potential problem with an expected increase in complexity. The equations governing 
the elastic displacement in the medium described above and depicted in Fig. [Tjis 

p.8) 

where body forces have been neglected. The associated Green’s tensor of the host phase 
satisfies 

as well as the far-field condition limx^oo Gij(x) = 0, noting that Gij = Gji. The resulting 
displacement field in the medium may be straightforwardly derive in integral equation form as 

f dG 

Ui{y) = u*{y) - {Clanki - clinu) “ y ) dx, ( 2 . 10 ) 

which holds for all y. Here u*(y) is the solution to the equivalent problem satisfying (12. 8p 
with no inhomogeneity present, or equivalently As in the potential problem, 

take derivatives of both sides of (|2.10p to form the strain tensor Cij = {dui/dxj + duj jdxi) j2, 
using the property dGki/dxj = —dGki/dyj so that we have, for all y £ V, 


('ijiy) ®ii(y) T {^mnk£ ^mnki) 




dyidy 


'j JVi 


^mn (x)Gfci(x-y) dx 


( 2 . 11 ) 




Here the notation indicates symmetry with respect to these indices, i.e. defining 


Qmnki — / 6mn(x)Gfcj(x y) dx. 


Wi 


we have 


d Qmnki 


dyedyj 




I/O Qmnki 0 Qmnli . 0 Qmnkj 0 Qmnij 
4 V OyiOyj dykOyj Oy^Oyi OykOyi 


( 2 . 12 ) 


(2.13) 


3 Uniformity of the Hill and Eshelby tensors 

3.1 The potential problem 

Impose so-called homogeneous temperature gradient conditions (in the language of microme¬ 
chanics, i.e. such conditions would lead to a homogeneous temperature gradient in an inhomo¬ 
geneous medium) in the far field, i.e. as |x| —>■ oo, 

T OiXi, (3.1) 

where 6i is uniform and therefore T* = OiXi and e* = 6i. Referring to (EID, one then asks, 
is there an inhomogeneity of any shape that can give rise to a uniform temperature gradient 
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field inside the inhomogeneity, i.e. for y G Vi? If such an inhomogeneity does exist, then (EZD 
is only consistent for y G Id if the tensor defined as 

is also uniform, i.e. is independent of y. The tensor P with components Pij defined in (13.2D 
is known as Hill’s Polarization (P) tensor for the potential problem and it possesses the 
symmetry Pij = Pji. If P is not uniform, it would mean that the assumption of a uniform 
temperature gradient field inside the inhomogeneity was incorrect. 

It transpires that when the inhomogeneity region is ellipsoidal the P-tensor defined in (13.2D 
is indeed uniform. This is proved in Appendix lA.ll where it is also shown that the general 
form for the P-tensor can be defined in terms of an integral over the surface of the unit sphere 

52. 


General form of Hill’s tensor for the potential problem: 
Ellipsoid in an unbounded medium 


The components of Hill’s tensor are defined as 


peiiipsoid_ det(a) 

ij - 4^ 


/ (f 13/2 


(3.3) 


where ^ is a unit vector that points from the origin, i.e. the 

centre of to its surface. Additionally is given by 










(3.4) 


and a is a second order tensor whose components are defined by 

3 

'*'■ '•'jn 


ttij — 'y ^ OndinSj 


(3.5) 


n=l 


as long as the semi-axes of the ellipsoid are aligned along the xi,X 2 and X 3 
axes, so that det(a) = 010203 . In fact it is always possible to define a in 
this manner by choosing xi,X 2 and X 3 to be aligned along the semi-axes 
of the ellipsoid, as long as one is happy for the principal axes of Cfj to be 
defined in different directions to xi,X 2 and X 3 should the principal axes of 
Cfj and Siij not be aligned. 


Clearly since the integral is over the surface of the unit sphere 5^, it is sensible to re¬ 
solve into spherical coordinates for the purposes of evaluating this integral. The form of 
(j3.3D illustrates the important general result that the P-tensor is uniform for an arbitrarily 
anisotropic ellipsoidal inhomogeneity embedded inside an arbitrarily anisotropic host phase. 
That Eshelby’s (weak) conjecture is true for anisotropic potential problems [36], [63|, means 
that the ellipsoid is the only shaped inhomogeneity for which the interior temperature gradient 
is uniform under all such far-held conditions of the form (13.ID . 

To determine the appropriate P-tensor in any circumstance then one can appeal to (13.3p 
and carry out the necessary integration. Alternatively, as shall be shown in 0 in many cases 
it is relatively straightforward to use symmetry arguments and results from potential theory in 
the isotropic host case together with scalings in some cases of host anisotropy, in order to derive 
explicit results, often in a more straightforward manner than directly evaluating the general 
result (|3.3p . In fact in the potential problem context, symmetry arguments and results from 
potential theory [38| are often sufficient to derive results for many special cases of ellipsoids 
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in host media that are at most orthotropic. The general result (I3.3p is thus suitable for more 
complex anisotropies than orthotropy or for example if the semi-axes of the ellipsoid are not 
aligned with the axes of symmetry of host anisotropy. 

We should recall that the P-tensor is independent of the anisotropy of the inhomogeneity 
and therefore we can retain arbitrary anisotropy for the inhomogeneity domain. The only 
aspects of the inhomogeneity that influence the P-tensor are its shape and, for anisotropic host 
phases, its orientation with respect to the axes of anisotropy of the host phase. It is important 
to note the following three points: 

• In the host region Vq the temperature gradient is generally not uniform. 

• For non-homogeneous temperature gradient conditions in the far field, the temperature 
gradient held inside an ellipsoidal inhomogeneity is generally not uniform. However if the 
prescribed temperature gradient is a polynomial of order n, then so is the held inside an 
ellipsoidal inhomogeneity, see [2]. This is known as the polynomial conservation property 
for ellipsoids. 

• Generally for non-ellipsoidal inhomogeneities in unbounded domains and general shaped 
inhomogeneities in bounded host domains V, the temperature gradient inside the inho¬ 
mogeneities is not uniform, although interacting E-inclusions |63] can lead to uniform 
interior strains and for specihc loadings, non-ellipsoidal inhomogeneities can yield uni¬ 
form interior strains, e.g. the counterexample of the Strong Eshelby conjecture given by 
Liu [63] . 

Regarding the hrst point, once we know the interior held ()3.16D . we can use this to determine 
the exterior held by using (12.61) so that for y ^ Pi, 

/dC 

^ ^(x-y) dx (3.6) 

where Aij is the temperature gradient concentration tensor linking the interior temperature 
gradient to that in the far-held, i.e. see T3.31 The gradient of (13.6p is not uniform since y 
now lies outside Vi. 


3.2 Elastostatics 


We used the symmetry relation Cijki = Cijik in deriving (12.lip as this turns out to be prefer¬ 
able in various contexts. Analogously to the potential problem, let us take homogeneous 
displacement gradient conditions in the far held, i.e. as |x| —oo 


Ui CijXj, (3.7) 

where Cij is uniform and therefore u* = CijXj and eL = (cjj -|- eji)/2. Note that eij does not 
have to be symmetric but if it is then it is simply the strain in the far held. As in the potential 
problem case the aim is then determine if there exists an inhomogeneity of any shape that 
is consistent with the assumption of uniform interior strain. If such an inhomogeneity exists, 
(12.lip is only consistent for y G Pi if the tensor dehned as 


Pijkeiy) 





(3.8) 


is uniform. The tensor dehned here is the P-tensor for elastostatics. It possesses the minor 
symmetries Pijke = Pijik = Pjikt by construction. Furthermore, thanks to the symmetry of 
the free space Green’s tensor Gij = Gji it also possesses the major symmetry Pijki = PkUj- 
It transpires that when the inhomogeneity region is ellipsoidal the P-tensor dehned in (13.81) 
is indeed uniform. This is proved in Appendix lA.ll where it is also shown that the general 
form for the P-tensor can be dehned in terms of an integral over the surface of the unit sphere 

52. 














General form of Hill’s tensor for linear elastostatics: 
Ellipsoid in an nnbonnded medinm 


The components of Hill’s tensor are defined as 


□ellipsoid_ det(a) 





^mn^np-^p 




dS{i) 


(3.9) 


where ^ and S'^ are as defined for the potential problem and a is defined 
in (13.5p . Furthermore 






where Nij is defined via 

N^kNuj = Nijil) = 


(3.10) 


(3.11) 


That Eshelby’s (weak) conjecture is true for isotropic elastostatics problems [36], [63] . 
means that the ellipsoid is the only shaped inhomogeneity for which the interior temperature 
gradient is uniform under all such far-held conditions of the form (|3.7I) . We stress however that 
it is not yet clear whether the weak conjecture is true in the context of anisotropic problems. 

To determine the P-tensor for an ellipsoid for a given host anisotropy one merely has to 
evaluate the surface integral in (|3.9p which can evaluated numerically very efficiently. For 
host anisotropies more complex than transversely isotropic it is generally recommended that 
the form (13.91) be employed and integrals are evaluated numerically. In what follows here the 
P-tensor shall be determined in the case of an isotropic host phase by appealing to various 
symmetries and potential theory. An important result derived by Withers m associated 
with a transversely isotropic host phase is also stated. 

As in the potential problem, the only aspects of the inhomogeneity that influence the P- 
tensor are its shape and, for anisotropic host phases, its orientation with respect to the axes of 
anisotropy of the host phase. Note also that the same three points described for the potential 
problem, preceding equation (|3.6p . also hold here in the elastostatics context. Furthermore, 
once the held is known inside the inhomogeneity region Vi the exterior field can be determined 
in terms of the Green’s tensor, as 

Ui{y) = ^ijVj - {Clankl - C^nkl)'^mnpqe*pg 

where Aijkt are the components of the strain concentration tensor (see ^3.3p . which links the 
strain inside the inhomogeneity to that in the far held. 

3.3 The Newtonian potential gradient and strain concentration tensors 

Defining the volume average 

7=1^^ /(x) dx (3.13) 

of the function /, it is straightforward to show that in the case of the conditions (|3.1I) . the 
body averaged temperature gradient is 


ei = ei = e*. (3.14) 

Of immediate interest is the temperature gradient field e\ inside an ellipsoidal inhomogeneity 
Vi, which from the theory developed above has been shown to be uniform so that it is equal 
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to its phase average, e\ = e\ where the phase average is defined as 

As such in the case of an isolated ellipsoidal inhomogeneities with homogeneous far-held con¬ 
ditions (j.S.ip . using (13.141) . the expression in (12.7p becomes 

e\ = e,-{Cl^-C%)e\p,,. (3.16) 

Using the symmetries C^- = Cji and Pjj = Pjj and re-arranging, (j3.16p can thus be written in 
the form 


= (3.17) 

Therefore one can relate the uniform temperature gradient inside the inhomogeneity to the 
average temperature gradient inside the entire body via a second order tensor, which is thus 
identihed as the temperature gradient concentration tensor for this problem. 


Temperature gradient concentration tensor: 

Ellipsoid in an nnbounded medium 

For an ellipsoidal inhomogeneity Vi embedded in an otherwise unbounded 
uniform medium, if homogeneous temperature gradient conditions (13.11) are 
prescribed in the far held, we have the exact relationship 

ej = el = AijCj (3.18) 

where the uniform concentration tensor Aij is dehned by 

AikAkj = dij, Aij = 6ij — C^j) (3.19) 

and Pij is dehned in (|3.2p . 


Note that Aij is the concentration tensor associated with an isolated inhomogeneity inside 
an unbounded host medium. The calligraphic notation Aij has been used to stress the link with 
(and distinguish from) the exact concentration tensor, usually dehned as Aij, and which links 
the phase average of the true temperature gradient inside an inhomogeneity to that in the far 
held in a complex inhomogeneous medium, which may consist of interacting inhomogeneities. 
For a dilute medium where interaction effects are not important, = Aij. 

Moving on to the elastostatics case, it is straightforward to show that in the case of the 
conditions O, the body averaged strain is 

Cij = P ~ ^ij- (3.20) 

The (uniform) strain ejj inside an ellipsoidal inhomogeneity Vi is thus equal to its phase 
average, ejj = e\y Therefore for an isolated ellipsoidal inhomogeneity with homogeneous 
far-held conditions (I32D, using (I3.20p . the expression in (12.lip can be used to determine the 
expression 


&ij — {lijkl + Pijmn{C 


(3.21) 


where we remind the reader that lijki is the fourth order identity tensor dehned in (IC.25h . 
Therefore the uniform strain inside the inhomogeneity can be related to the average strain 
inside the entire body via a fourth order tensor, which is thus identihed as the strain concen¬ 
tration tensor for this problem. 
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Strain concentration tensor: 
Ellipsoid in an nnbonnded medinm 


For an ellipsoidal inhomogeneity Vi embedded in an otherwise unbounded 
medium, if homogeneous displacement conditions (|3.7p are prescribed in 
the far field, we have the relationship 


(3.22) 

where the uniform concentration tensor Aijke. is defined by 

•^ijmri'^mnki — lijkit — ^ijki ^mnki)' (3.23) 


4 The potential problem: specific cases 


4.1 Isotropic host phase 


Assume that the host phase is isotropic, so that C?- = ko6ij and therefore the associated 
free-space Green’s function is 


From (13.2p and (|l.ll) therefore 


G(x-y) 


1 1 

dvrfeo |x - y 




1 

ko dxidxj ’ 


S'jj(x) 


d^F 

dxidxj 


where F is the potential defined by 


r(x) 


An Jy^ |x - y 


dy. 


(4.1) 


(4.2) 


(4.3) 


Note that this is the negative of the Newtonian potential (see for example Kellogg [48]) asso¬ 
ciated with an ellipsoidal domain Vi. From potential theory therefore 


v2r(x) 


d^r 

dxkdxk 



(4.4) 


and furthermore F(x) is a quadratic function of the components of x (see Appendix [B|), illus¬ 
trating the uniformity of the P-tensor in this case. 

As an aside, note that since the host phase is isotropic, the temperature field exterior to 
the inhomogeneity is determined via (13.6p . i.e. 


Tiy) = OiVi + (C^. - ko5kj)Akee£^^^^. (4.5) 

This solution tends to Oiin in the far field, as it should do. 

Once Pij is determined for an isotropic host phase the associated concentration tensor for 
an isolated inhomogeneity may then be found from (|3.19D as is now illustrated in a number of 
special cases of specific inhomogeneities with given shape and anisotropy. 


4.1.1 Sphere in an isotropic host phase 


When Vi is a sphere, it is clear from (|4.3I) that F(x) must be spherically symmetric and hence 
. must be isotropic (and uniform), i.e. 


92 F 

dxidxj 




(4.6) 
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for some constant 7 . Note that the form ()4.6|) is a result of the spherical shape and not any 
assumption regarding isotropy of the inhomogeneity as such an assumption has not been made. 
Performing a contraction in (|4.6I) and using (|4.4p with x G Vi yields 7 = |. Therefore from 
(lOI) and (fTTH 


For practical purposes, especially for use in micromechanical methods for bounds and estimates 
of effective material properties, it is useful to write down the associated concentration tensors. 


Isotropic sphere 

If the spherical inhomogeneity is isotropic with conductivity Cb = ki5ij, one can show straight¬ 
forwardly using (13.191) and properties of second order tensors (see Appendix 1C.3p that 


jAi 


3fco 

ki + 2ko ' 


(4.8) 


Anisotropic sphere 

Consider a transversely isotropic sphere where the plane of isotropy is the X 1 X 2 plane. The 
conductivity tensor therefore takes the form = ki{Qij -|- v6i^6js) where Qij is defined 
according to 


^ij — ©ijr T 


(4.9) 


and V indicates the degree of anisotropy, with v = 1 giving isotropy. Using the result derived 
in (BZI) and (BJl) together with properties from Appendix 1C.31 the concentration tensor Aij 
can be written down in the form 


_ 3fco 
ki + 2kQ 


Qij 


3ko 


vki + 2ko 




(4.10) 


Setting V = 1 recovers the isotropic result (14.80 . 

Averaging over all orientations of the anisotropy of the inhomogeneity will yield an isotropic 
concentration tensor of the form 




(4.11) 


where the underline denotes averaging over orientations. By performing this orientation aver¬ 
aging (see Appendix 1C. 3. 4p on (I4.10p it is straightforwardly shown that 


2kQ 

k\ -\- 2fco 


kp 

vki + 2kp 


(4.12) 


4.1.2 Circular cylinder in an isotropic host phase 

When Vi is a circular cylinder with axis of symmetry in the X 3 direction, it is clear that r(x) 
should be independent of xp and isotropic in the X 1 X 2 plane so that 


d^T 

dxidxj 


70b 


(4.13) 


for some constant 7 where Qij was defined in (14.91) . Performing a contraction in (I4.13p and 
using (|4.4p the result 7 = ^ is obtained. From (|4.2p therefore 


p.. —_P) .. 

“ 2 fco 


Sij = 



(4.14) 
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If the cylinder is isotropic with conductivity tensor 


Clj = 


ij, 


(4.15) 


it is straightforward to show that 

•^ij = ^ ^ ©ij + (4.16) 

Alternatively, suppose that the cylinder is transversely isotropic with conductivity tensor 


Clj = h{Qij + vdisSjs). 


(4.17) 


Interestingly one can show that in this case the concentration tensor is identical to the isotropic 
case, i.e. that in (I4.16p : the parameter v does not appear in the concentration tensor. Of course 
if the axis of symmetry of transverse isotropy is not aligned with the cylinder axis then this 
concentration tensor would then depend on v. 

If the (uniform) orientation average of (|4.16p is taken, the associated concentration tensor 
is derived: 


_ ki+5ko 

-*^'“3(fei + feo) 


(4.18) 


This last result is often used when very long, thin needle-like inhomogeneities are uniformly 
distributed and oriented throughout some host medium. 


4.1.3 Ellipsoid in an isotropic host phase 

Consider now the general case of an ellipsoidal inhomogeneity and as usual denote the semi¬ 
axes of the ellipsoid as aj, j = 1,2,3. It is straightforward to show, using the theory of the 
potential, as in Appendix iBl that for an ellipsoid in an isotropic host phase, the function r(x) 
is quadratic in the components of x and can be written in the closed form 


/X? Xn Xo \ dT 

r(x) =4 + ^ + ^- it-V 

V af ai J ^ aj daj 


where 


T = 


1 r 

-010203 

In Appendix [B] it is then shown that 


dt 


y/{a( + t){a^ + t){al + t) 


dxidxi — 

n=l 




where with = 03/an-, 

roo 

£{x;£i,£2) = Y Jq 


ds 


(1 -h sx2)^(l -t- sef)(l -h s£\){l -h s) 


Therefore 


(4.19) 


(4.20) 


(4.21) 


(4.22) 


P — 


1 

fco 


3 

n=l 


Sij — k()Pi 


ij- 


Finally note that using (|4.4h it is easily shown that 71 -|- 72 + 73 = 1- 


(4.23) 
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4.1.4 Spheroid in an isotropic host phase 


Denote the semi-axes of the spheroid as ai = 02 = a 7 ^ 03 and use this in ()B.27p which becomes 


T = -a-. 


"L 


cos V’ 


0 -I- (1 — £2) -0 


(I'll! 


where e = 03 /a. Make the substitution (3 = sin^ to find 


T = -ai 


dp 


£2 + (1 - £ 2)^2 




arccosh(e:) 

1 ’ 

arccos(£’) 

£ > 1 

£ < 1 

1, 

£ = 1 


noting that £ = 1 is the case of a sphere. 
Therefore from (I4.19p it is clear that 

1 


r(x) = ^ixi + xl)T{e) + -xlS{e) - T 


(4.24) 


(4.25) 


(4.26) 


where 


5(e) = 4 t- 

“3 


2 dT 
03 dos 


1 £ 

1 — £2 1 — £2 


. Vl— 


arccosh(£), 

arccos(£), 


£ > 1 , 

£ < 1 


(4.27) 


and T(£) = ^(1 — 5(£)). The function 5(£) has taken many forms in the literature but it is 
felt that this is a most clear, consistent and concise formulation. Note that 5(£) —^ as £ ^ 1 
for the spherical case (see further details below). 

Using (14.2611 in ()4.2I) . the resulting Hill and Eshelby tensors take the form 


Pij — ^ T ) 


Sij = koSij (4.28) 


where 73 = S{e) and 7 = T(£). 

Expressions for the concentration tensors associated with the spheroidal inhomogeneity 
case can now be determined straightforwardly. Eor an isotropic spheroid, 

ko + [ki - koh ko + {ki - kohs 
Averaging uniformly over orientations of the axes of the spheroid yields 


A-■ 

—ij 


1 / __feo_\ ^ 

3 \ko + {ki- koh ^0 + (^1 - ^ 0)73 ) 


(4.30) 


One can take limits in the case of the spheroidal inhomogeneity in order to derive the 
following results, some of which confirm cases considered above. When Vi is 

(i) a sphere, i.e. £ ^ 1 , it is deduced that 7 = 73 = |, 

(ii) a cylinder, i.e. e —>■ 00 , it is deduced that 73 = 0 , 7 = ^, 

(iii) a disk or layer, i.e. £ —)• 0 , it is deduced that 73 = 1 , 7 = 0 . 


When used in (|4.29|1 (i) and (ii) confirm the results derived for the concentration tensors for 
isotropic spheres and cylinders derived in S ^4.1.1I and 14.1.21 respectively. One has to be rather 
careful in taking these limits and for (i) use L’Hopital’s rule appropriately, noting that as £ ^ 1 

5(£) = - - —(£ - 1) + —(£ - 1)^ -I- 0((£ — 1)^). (4.31) 

In (ii) one has to use the fact that arccoshx ~ logx as x ^ 00 . The result for layers in (iii) 
can also be obtained via straightforward symmetry arguments. 
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4.1.5 Limiting case of an elliptical cylinder 

One can use the formulation for general ellipsoids above in order to obtain a result for an 
elliptical cylinder, unbounded in the X 3 direction with semi-axes ai and 02 lying along the xi 
and X 2 axes respectively. Taking the limit 03 —)• 00 in ()B.32p . one can show that 


02 r 

dxidxj 


^ ^ 'yn^in^jn 


n=l 


where 


71 = 


ai 02 


/*CXD 

Jo (aj 


ds 


+ s)^{al + s )2 


72 = 


0102 


ds 


'0 (of -I- 5)2 (af -I- 5)2 


(4.32) 


(4.33) 


noting the fact that X 3 dependence is eliminated as should be expected. The integrals can be 
determined explicitly, noting the indefinite forms 



(of-|-s) 2 (a|-|-s) 2 (®2 ®l) + ® 


of -1- S V 

of + s j ’ 

(4.34) 

of -b S \ 2 

02 + sj 

(4.35) 


and therefore 


71 = 


02 


oi + 02 

The Hill and Eshelby tensors therefore take the form 

1 


72 = 


Ol 


Ol + O 2 


p = - 

ko{ai+a 2 ) 


(o2(5ji(5ji -|- aidi 26 j 2 ), 


Sij — koPi^j 


(4.36) 


(4.37) 


As regards the concentration tensor for an isotropic cylinder with Ck = ki5ij, this is determined 


in the form 


^0(1+ e);; ;; , fco(l + e);; ;; , X X 

;- OilOjl -|- — -—- di2dj2 + O^dj^ 


ko + kie ' ki + koe 

where e = 02/01 is the aspect ratio of the ellipse. 


(4.38) 


4.1.6 Limiting cases of a cavity, penny-shaped crack and ribbon-crack 

It does not really make sense to define a temperature gradient concentration tensor in the 
context of cracks or cavities because clearly there is no interior field. However it turns out that 
this concept is useful and can be interpreted as linking the far-held to the held on the surface 
of such inhomogeneities m with an appropriate dehnition of “cavity temperature gradient”. 
As such here the results above are used in order to derive associated concentration tensors for 
cracks and cavities. 

Consider a spheroidal inhomogeneity and the limit /si —)• 0 in (I4.29p . This yields 

jAij = - &ij + Tj (4.39) 

1-7 1-73 

This is the concentration tensor for potential problems involving spheroidal cavities. 

Next consider the so-called “penny-shaped crack” limit. We require the asymptotic form 
of 7 (e) and 73 (e) as e ^ 0. These are easily shown to be 

73 (e) = 1- |e + 2e^ + 0(e^), 7 (e) = -F O(e^). (4.40) 
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As such one can derive the form 

•^ij = Qij + 2 ^ ^i 3 ^j 3 + 0{e) (4-41) 

where expansions have been taken for e <C 1 and terms up to 0 ( 1 ) have been retained since 
higher order terms will clearly vanish as e —>• 0 . 

The coefficient of 6136 in (I4.4ip involves an apparently singular limit as e —>• 0. That 
this is not a problem arises from the fact that this expression is used in formulae for effective 
properties of cracked media where this term is always multiplied by a volume-fraction term 
in such micromechanical methods, (or rather a “crack-density”) that is proportional to £ [38], 
[39] . Note that taking the limits in the opposite order, i.e. e —)■ 0 and then fei —)• 0 yields an 
inconsistent result, giving rise to singular effective material behaviour in the crack limit which 
cannot be correct. 

Finally, consider a different limit, the so-called “ribbon-crack limit”. Take ki = 0 in the 
elliptical cylinder result (|4.38p to find 

= (1 + + hishjs- (4.42) 

Therefore as e ^ 0 

A-ij = —Si 26 j 2 -l- Sij 0(e). (4.43) 

As in the penny-shaped crack result above, the concentration tensor for the ribbon-crack is 
singular. 


4.2 Anisotropic host phase 

The general form (13.3p for the P-tensor associated with arbitrary host anisotropy requires the 
necessary surface integral to be evaluated. In the case of transversely isotropic and orthotropic 
media however, where principal axes are aligned with the semi-axes of the ellipsoid, the problem 
can be simplified significantly by employing a scaling of the Cartesian variables in order to 
reduce the isolated ellipsoidal inhomogeneity problem in an anisotropic medium to the case of 
an ellipsoidal inhomogeneity (with different semi-axes) in an isotropic medium. Therefore the 
results derived above for the isotropic host phase case can be used in the scaled domain and 
then map back to the physical domain to obtain the appropriate physical Hill and Eshelby 
tensors. 

As usual consider the case of an ellipsoid with semi-axes aj,j = 1,2,3 but now embedded 
in an orthotropic host medium (with principal axes aligned along Xj, i.e. with the semi-axes of 
the ellipsoid) so that 

Oij = ko + V 26 i 26 j 2 -I- vsdisdjs ), (4.44) 

where U 2 = 1 (or U 3 = 1) for transverse isotropy. The governing partial differential equation is 

A A A 

dxi \ dxj 

Now employ the simple rescaling 



j = 1, 2,3 (4.46) 

where ui = 1 is introduced for notational convenience (the conductivity along the xi axis is thus 
ko), so that the semi-axes of the ellipsoid in the mapped domain become aj = aj / y/v], j = 1,2,3 
and denote the scaled ellipsoid as Vi. The governing equation then becomes that governing 
isotropic media so that 




a^r(x) 

dxi 


^in^jn — 


1 ^ 

-A 

ko V„ 


1 a2r(k) 


dxi 


^inhjn 


(4.47) 
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where r(x) is defined in terms of the isotropic (due to scaling) Green’s tensor as defined in 
(14.ip but now integrated over the scaled ellipsoid Vi, i.e. 


f(i) 


1 

dvr 




(4.48) 


As a consequence for a general ellipsoid, the result (|B.29p can be used but with Xj replaced by 
Xj and Oj replaced by aj,j = 1,2,3. Therefore with reference to (|4.22l) 


d^T ^ ^ 

dx dx ~ (4.49) 


where in = cl's/an- The P and S-tensors for an ellipsoid embedded inside an orthotropic host 
medium can then be written as 


P — 




in^jni 


n=l 


3 

Sij — ^ ^ fJn'yn^in^jn (4.50) 

n=l 


where 7 ^ = (e^;ei, 62 ). 

Note that the above scaling approach is considerably simpler than carrying out the nec¬ 
essary integrals in the corresponding general expression (13.3p for the P-tensor. Now consider 
some specific cases of anisotropy of the host phase. First consider an inhomogeneity embedded 
in a transversely isotropic (TI) host phase with conductivity tensor 

Cfj = koiQij -I- vdisdjs). (4.51) 

Hence the (orthotropic) P-tensor for an ellipsoid with semi-axes aligned with the principal 
directions of anisotropy is found by setting vi = V 2 = 1 and U 3 = u in (|4.50l) above. Simplifi¬ 
cations arise for a spheroid of course as shall now be illustrated. 


4.2.1 Spheroid in a transversely isotropic host phase 

Consider a spheroid in a transversely isotropic medium where the major/minor axis of the 
spheroid is aligned with the axis of transverse isotropy of the host phase. Denote the semi¬ 
axes of the spheroid as a = oi = 02 / 03 and the axis of transverse isotropy as X 3 . We use the 
scaling argument above to see immediately that the P and S-tensors are given by 

Pij = ^ {jQij + , Sij = koPij (4.52) 

/Co 

where with reference to (I4.27p 

73 (e) = ^5 , (4.53) 


e = 03/0 and 7 = ^(1 — U 73 ), the latter being derived by using 


d2f Q2p Q2p 

dxf 8x2 dxl 

for X G Pi. As is evident, given the calculations already made associated with isotropy, this 
is a much simpler mechanism for obtaining results for an anisotropic medium than using the 
general form for the P-tensor and carrying out the necessary subsequent surface integral. 

Assuming the spheroid itself is isotropic with conductivity tensor = ki6ij, and using 
(I4.5ip together with the form of P-tensor in (|4.52p . the concentration tensor defined in (I3.19P 
can be straightforwardly determined as 



3..- feo 

ko + {ki - koh 


Qij + 


ko 


ko + {ki - vkohs 




(4.55) 
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with 7 and 73 as defined above. Alternatively, supposing that the spheroid is now transversely 
isotropic with the same axis of symmetry as the host, i.e. Cb = ki{Sij + one finds 

that 


A.. = _^ 

ko + {ki - ko)'y 


Qij + 


kn 


ko + {(ki - vko)j3 




(4.56) 


4.2.2 Circular cylinder in a transversely isotropic host phase 

The circular cylinder limit can be taken in the spheroid case considered in 114.2.11 where the 
cross-section of the cylinder sits in the plane of isotropy of the host medium. It is then 
anticipated that the P and S-tensors will be TI. It has been discussed above that 5(x) —)• 0 
as X —)• 0 and therefore as with the isotropic host case from (j4.53D 73 —>• 0 . As such 7 = 
^(1 — 1 ^ 73 ) = ^ and then 

^ij — ~ 2 ^*-^ 

so that in fact this tensor is unchanged from the case of an isotropic host phase as in (j4.14D . 
The concentration tensor for an isotropic cylinder can be straightforwardly determined as 

•^ij = ^ ^ ©ij + ^i3^j3- (4.58) 

The concentration tensor associated with a transversely isotropic cylinder is also given by that 
in (14381) . 

An interesting non-standard example is the case of a spheroid embedded inside a trans¬ 
versely isotropic host phase where the axes of symmetry and semi-axes are not coincident. In 
this case the general (surface integral) form of the P and S-tensors must be used with the semi¬ 
axes aligned with the x axes but with all components of the modulus tensor being generally 
non-zero. 


4.2.3 Ellipsoid in an orthotropic host phase 

Consider an ellipsoid with semi-axes aj,j = 1,2, 3 that are aligned with the axes of anisotropy 
of the host medium with orthotropic conductivity tensor as defined in (j4.44p . Analogous scaling 
arguments can be used as above in order to scale this problem into an ellipsoid in an isotropic 
host and then scale back to the physical domain, as described above to show that 


1 ^ 

Pij — ^ ^ 


jn 


(4.59) 


n=l 


n=l 


where 


dn — ^ {^n ) £1) £2 ) 


Vr, 


(4.60) 


where ui = 1 and noting that in = where 

For reference, P-tensors for a variety of problems are summarized in table [TJ 
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Host anisotropy 


Isotropic 



Transversely 

isotropic 

Vl = V2 = I V'i = V 



Orthotropic 

Vi = l^V2i^V3 


Worse than orthotropic 

or semi-axes of ellipsoids 

not aligned with axes 
of anisotropy. 


Inclusion shape 


Ellipsoid 
«i / 02 7^ 03 
£fi — O3/ Chn 


Spheroid 
oi = 02 = a / 03 

e = a^/a 


Sphere 


Ellipsoid 
Or / 02 7^ 03 

£n = 03/af), 


Spheroid, ai = 02 = a / 03 
and 03 is aligned with 

axis X 3 of transverse isotropy 


Spheroid, ai = 02 = a / 03 
and a is aligned with 
axis X 3 of transverse isotropy 


Sphere 


Ellipsoid 
Or 7^ 02 7^ 03 


P-tensor 


Use potential theory; 

1 3 

Pij — ■; ^1 1 

fco 


Use potential theory; 
Pij — j T 

7 = x(l - 73 ), 73 = 5(e) 


Use symmetry; 

P • — d- • 
“ 3A;n 


Use scalings and potential theory; 
Pij = X]n=i :^5(en; ei, e2)<Jm<Jjn 
— 03 /h^ and Qjyi — Ojrfi!PVri‘ 



Use scalings and potential theory; 

1 3 1 . . . 

Pij = -r- X]n=i :^5(en; ei, e2)<Jm<Jjn 

ft-O 

— 03 / 0 . 77 , and O 77 , — ^nl\fPn' 




Use scalings and potential theory; 

Pij = Z]n=l :^5(en; ei, e2)<Jm<Jjn 

— O 3 /O 77 , and O 77 , — ^n!\fPn‘ 


Use general integral form; 
^ _ det(a) r ^ij dS(i) 
dvr Js '2 

^^J=U,/{C^kMe) 


Table 1; Table summarizing results for the P-tensor associated with ellipsoidal inhomogeneities 
for potential problems. 
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5 Elastostatics: specific cases 

5.1 Isotropic host phase 

For the case of an isotropic host phase case the elastic modulus tensor is defined as 


^ijki 


(5.1) 


in terms of the isotropic fourth order basis tensors (IC.28h and (IC.24h . Here kq and fiQ are the 
bulk and shear moduli of the host, noting the relation to Poisson’s ratio uq 


_ 2 /ro(l + i^o) 
“ 3(1 - 2uo) ■ 


(5.2) 


The appropriate isotropic Green’s tensor is 


Gy(x-y) = 


1 


1 


52|x - y| 


Att/jq |x — y| 167r//o(l — z^o) dxidxj 
and the expression for the P-tensor in (|3.8I) therefore becomes 


(5.3) 


„ 1 f d^r ^ d‘^r ^ d‘^r ^ \ 

^ijki — - I 7^ 7 —Oik + "7 7 On + — —Ojk + -7 7 Ojl 1 

4/io \OXjOX£ OXjOXk OXiOXi OXiOXk / 

1 5% 

4^0(1 ~ ^o) dxidxjdxkdxi 

The potential P is that already encountered and defined in (|4.3I) and T is defined by 


(5.4) 




-y| dy, 


which satisfies (see for example Kellogg |48] 1 

V% = -2V'^T = - 2 x^(x) 


(5.5) 


(5.6) 


Once Pijki is determined, the components of the Eshelby tensor can be calculated from (II.ip 
and the associated concentration tensor can be found from (|3.23p . Recall that no assumptions 
have been made regarding the anisotropy of the inhomogeneity. This is not required in order 
for the P-tensor to be determined. The only aspects of the inhomogeneity that influence the 
P-tensor are its shape and, for anisotropic host phases, its orientation with respect to the axes 
of anisotropy of the host phase. 


5.1.1 Sphere in an isotropic host phase 

Assume that the host phase is isotropic with elastic modulus tensor given in ()5.ip and consider 
the case where Vi is a sphere. The (uniform) tensors d‘^T/dxidxj and /dxidxjdxkdx^ will 
be spherically symmetric, i.e. isotropic and must possess full symmetry with respect to the 
interchange of any index. As with the potential problem therefore, 


(9^r _ 1 
dxidxj 3 


(5.7) 


and for T the general isotropic form, with the additional constraint regarding symmetry with 
respect to interchange of any indices, is 


"3 3 3 3 '^{.dijdki T hikdji Siidjk') 

OXi ox j ox k OXfi 


(5.8) 
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where V’ is a constant to be determined. Performing the contractions j = i, (. = k in (15.8p and 
utilizing ()5.6I) . V' = —2/15. Using (|5.7p and (15.8p in ()5.4I1 . the components of the P-tensor are 


“ 6//0 ^ 30/io(l - yo) 




(5.9) 


After simplification and writing in terms of the tensors becomes 

Pijkt = Pilfjke + P^lfjke^ ( 5 - 10 ) 

where upon using the relationship (15.21) the components of the P-tensor are determined as 


_ 1 - 2z/o _ 1 _ 4 - 5zzo _ 3(Ato + ‘^Po) 

6po{l - vq) 3ko + 4/io ’ 15/xo(l - i^o) 5^o{3ko + 4|Uo) ’ 

From this form the Eshelby tensor is easily determined via dni) as 

Sijki = silijki + S2lijkej 

where 

_ 1 + zzq _ 3ko _ 2(4 - 5j/o) _ 6(ko + 2/Uo) 

3(1 - zzo) 3Ko + 4/ro’ 15(1-1/0) 5(3 ko+4/xo)‘ 


(5.11) 


(5.12) 


(5.13) 


Isotropic sphere 

The P-tensor derived in (|5.10p holds for a spherical inhomogeneity of arbitrary anisotropy, 
embedded inside an isotropic host phase. If the inhomogeneity is also assumed isotropic with 
elastic modulus tensor of the form 


'^ijki — + ‘^Pl^ijki^ (5-14) 

one can use ()C.25p and the contraction properties of the tensors lljki, ifjke defined in Ap¬ 
pendix [U31 in order to write 


■^ijke — hjki 4“ hjkl + PijrriniCmnkl ^rnnki) 

= (1 3(ki - Ko)pi) lljki + (1 + 2(Aii - Po)P2) lijkt (5-15) 

Finally the inversion properties of such a tensor are employed (as described in Appendix [CAiS]) 
together with (15.111) in order to obtain the strain concentration tensor: 

_/3ko_+^\ _ 5/io(3Ko -H4^o) _ ^ 

V3«i+4//oy 3^o(3/io + 2m)+4/.o(2/io + 3w) *4^^- ^ 

Taking ki,/xi —>• 0 in (I5.16p the concentration tensor for a spherical cavity is obtained purely 
in terms of the host Poisson ratio, i.e. 


•^ijkt — 


3(1-^o) , 

2(1 - 2 i/o) 


15(1 - yp) j2 

7-51/0 


(5.17) 


Cubic sphere 

Assume now that the sphere has cubic symmetry with elastic modulus tensor 

Cijki = ^'^l^ijki + ‘^Pl^ijke + Vl^ijkE 

where 


^ijkl 



i = j = k = e, 
otherwise. 


(5.18) 


(5.19) 
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Instead of (I5.15j) . it is found that 


Ajki = (1 + 3(ki - Ko)pi) l}jkl + (1 + 2 (/Xl - Po)P2) ifjki + rnPijmnSmnkl- (5-20) 

Next, using the form of the P-tensor in (15.101) and the relationships (|C.85p . the expression 
(j5.2UI) becomes 

AjH = + ^^ifjki + ^3^ijk£, (5-21) 

where 

di = 1 + 3 (ki - Ko)Pi + ??i(Pi-P 2 ), q ;2 = 1 + 2(/ii -/xo)p 2 , ds = W 2 - (5.22) 

Finally, appealing to the theory regarding cubic tensors in Appendix IC.4.21 the concentration 
tensor is 


■^ijke — Oillijki + ^^ifjki + OiS^ijke 


(5.23) 


where 


ai = 


02 + Olds 


0^2 — 

02 


03 

0^3 . 

On 


02(di 03) Lx'z ^2 

Note that when r/i = 0 the case of an isotropic sphere is recovered as in (|5.15l) . 


(5.24) 


5.1.2 Circular cylinder in an isotropic host phase 

Now assume that Vi is a circular cylinder with axis of symmetry in the X 3 direction. As for the 
potential problem, r(x) must be independent of X3 and isotropic in the X 1 X 2 plane. Hence, as 
described in (I4.13p 


d'^T 

dxidxj 



(5.25) 


Furthermore d' must also be independent of X 3 and be isotropic in the X 1 X 2 plane and hence 
write 


d'if 

dxidxjdxkdxi 


^(.^ij^ki 4 “ ®ik®ji 4 “ ®il®jk\ 


(5.26) 


for some constant ip, where the fact that this tensor should be fully symmetric with respect to 
interchange of any of its indices has been used. Performing the contractions j = i, i = k and 
using (|5.6p leads to the conclusion that ip = —j. From (|5.4I) therefore 


Pijki — o i®j£^ik 4“ Qjk^i£ 4“ Qi£5jk 4“ Qik^j£^ 

OPO 

Y 6 /i~(l 4“ Qik®j£ 4“ Qi£Qjk)- (5.27) 

By writing dij = Qij + 6136 in the first term of (I5.27p and then recognizing the appropriate 
Hill TI basis tensors (see Appendix IC.4.3h that arise as a result of the various contraction 
terms, the P-tensor can be written in the form 


Pijki 


8po 


1 

4/io 


2(1 - 2120 ) 

I -1^0 


P'ijk£ 


+ 


3 - 4t/o 
I-uq 


'1/5 

^ijki 


2 l^o Ao + 3^ 5 

Ao + 2/xo Ao + 2^0 


4“ ‘^'P-ijki 
4“ 'H%k£ 


(5.28) 

(5.29) 


In order to derive the Eshelby tensor one contracts P with the host modulus tensor C^. 
In order to do this write the isotropic basis tensors in terms of the Hill basis tensors using 
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the expressions (IC.42p - ()C.44p and then use the contractions summarized in Table [2] in order 
to determine that 


Sijki — + (3 — + 2(1 — i^o)'H^jk^ ■ (5.30) 

It may well be the case that a different basis set should be used if it transpires that the cylinder 
itself is more anisotropic than TI, for use in concentration tensors for example. However for 
computation, the matrix formulation, as described in Appendix 1C.51 can be of great utility 
when anisotropic basis tensors are becoming rather cumbersome. 

Suppose now that the cylinder is isotropic with elastic modulus tensor as defined in (15.1411 . 
Once again using (|C.42p - (]C.44l) . constructing is then just a matter of exploiting the 

contractions in Table [2j It transpires that Aijki takes the form 


where 


+ /To 

(Ao + 2/io) 


OtA = 1 , 


6 

Aijki. = ^n'Hijf.i 

n=l 


(5.31) 


(Ai — Aq) 

2(Ao + 2^o)' 

1^0 + /ii(3 - 4t/o) 

4//o(l - i^o) 


0.3 = 0 , 


(5.32) 


de 


(m + Mo) 

2^0 


(5.33) 


Using the inversion expression for transversely isotropic tensors stated in Appendix 1C. 4. 31 the 
appropriate concentration tensor is thus derived as 


6 

= (5.34) 

72=1 

where 


Ao + 2/io 

(Ao — Ai) 

03 = 0, 

(5.35) 

(Ai + fii + fio) ’ 

2(Ai + + fio) ’ 

Oa = 1, 

_ 4/ro(l - 1 ^ 0 ) 

IJ .0 + fii{3 - Avo)’ 

2^0 

«6 — , 

+ Mo 

(5.36) 


Since 02 A “3 if is seen that ^1133 AI 3311 . The circular cylindrical cavity limit can be 

obtained by setting //i = 0. Using the expression A = 2^i//(l — 21/), 


2(1 - Vo) 
(1 - 2vo) 
Oa = 1 , 


«2 = «3 = 0, (5.37) 

1 — 2Vo 

as = 4(1 - vq), oe = 2. (5.38) 


Finally note that the fourth order tensor (15.3411 can be represented in matrix form (see Ap¬ 
pendix [C3]) as 


[Al] = 


where 


an — 


(1 - t^o)(3 - 4:Vo) 
(1 - 2vo) '■ 


/ 

On 

012 

013 

0 

0 

0 

\ 


012 

On 

013 

0 

0 

0 



031 

031 

O33 

0 

0 

0 



0 

0 

0 

033 

0 

0 



0 

0 

0 

0 

033 

0 


V 

0 

0 

0 

0 

0 

066 

/ 


012 = 


{vq - 1)(1 - 4 t / o ) 

(1 - 21^0) ^ 


ai3 = 




I- 2 A 

fl66 = 2(1 — Vq). 


(5.39) 


031 — 0 , 


033 = 1 
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(5.40) 

(5.41) 



































Clearly it is possible to write down explicit expressions for the concentration tensor when 
the circular cylinder is anisotropic. This merely complicates the tensorial (or matrix) oper¬ 
ations after the derivation of the P-tensor in (I5.27p . Given this P-tensor, perhaps the most 
important aspect is then to choose the tensor basis set correctly, given the anisotropy of the 
inhomogeneity. For example, when the cylinder itself is transversely isotropic (a common oc¬ 
currence in applications) it is considered sensible to use a TI tensor basis set. For practical 
purposes and especially for the sake of computation, using the matrix formulation of tensors is 
advantageous in cases where the tensor basis sets become rather cumbersome. The following 
procedure is used, in the usual notation, referring to Appendix 1C.51 and defining the 6 x 6 
matrix [P] associated with the P-tensor, define 

[A = [I] + [P][W][C^-C^] (5.42) 

where [W] is defined in (|C.57p and therefore 

[A] = [W]-^[A]-^[W]-\ (5.43) 


5.1.3 Spheroid in an isotropic host phase 

When Vi is a spheroid, the potential theory outlined in Appendix [B] is once again of use. It is 
clear that the P-tensor must be transversely isotropic and therefore will take the form 

6 

Pijke = (5.44) 

n=l 


The separate contributions to the P-tensor shall therefore first be decomposed into this form. 
Firstly, from the potential case described in Example 14.1.41 


92 F 

dxidxj 


'l&ij + IsSisdjs 


(5.45) 


where 73 = S{e) and 7 = ^( 1 — 73 ). Representing all terms in the Hill basis one can show that 


4 \dxjdx£ dxjdxk dxidxe ^ dxidx^ 

= liP-ijke + P^ki) + IsP-ijke + 2 (^"1 '^‘i)'P%kt 

Sik — 5 i 3 (ffc 3 and writing for example, 
d^F 

o o -(©ifc + 

C/ JL j (J JL £ 

{iQjl + 73<^i3<5£3)(0ifc + Sisdks) 

'fiQjeQik + QjAsha) + isi^ik^js^es + ^is^js^ka^es)- (5.46) 

Doing this for each term on the left hand side of ()5.46p . combining and using the definitions 
of the TI basis tensors in App endix I C.4.3 1 leads to the form on the right hand side of ()5.46p . 

Further, after much algebraic manipulation using the simplifications of the integrals in 
Appendix [B] in the case of spheroids one can show that 


This is seen by using = 
d^F 


dxjdxi 


^ik 


1 5% 

4 dxidxjdxkdxi 


6 

i’nTPijkl 

n=l 


where, upon using 73 = 1 — 27 , 


V’l 

IpA 


£ 2 ( 47 - 1 ) -7 
4(1-£2) 

37-1 

2 ( 1 ^’ 


'4>2 = V’3 



£ 2(1 - 27 ) -7 
4(1 -£2) 


(5.47) 


(5.48) 

V'e = 2V’2. (5.49) 
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Therefore 

pi = — \ i + Tr~—’ 

/ro V (1 - ) 

P 4 = — ( 1 - 27 + -- rlpi ] , 

Po \ (1 - z^o) / 

A good check is to ascertain that the result for a sphere in an isotropic host phase is recovered 
by taking e ^ 1 and using ()4.31l) . This is easily done and yields the result derived in T^.l.ll 
associated with a sphere. Another useful limit is to take e ^ 0. Results already derived can 
be employed, e.g. (j4.40p which when used in (|5.48p - (j5.49h yield 

2 

+ 0{e^), ^2 = zAs =+^+ 0(e^), (5.53) 

lb lb 2 


P 2 =P 3 = 


11^2 


/Uo(l - 1 ^ 0 )’ 

P 5 = — [1 + 7 T 71 -rVl 

Po V 2(1 - vq) 


(5.50) 

(5.51) 

(5.52) 


+—e — 2£^ + 0{e‘^), '4’b —+ 0{£^)^ ———e + + 0{£^). (5.54) 

For the components of the P-tensor this then gives 

(5.55) 

(5.56) 

(5.57) 

(5.58) 

(5.59) 

The Eshelby tensor with respect to a TI basis, in the form 


Pi = 

Pb = 

P& = 


Po 

U6(l-z.o)^ 

1 1 

( 

Po 

y 16 ( 1 -z/o) 

1 1 

( 1-2120 

Po 

^2(1-z/o) 

Ti 

( zi '(7 - 8120) _ 


£ + 


+ 0{e^), 

1 


+ 0 (£^), 


-£ — 


£M+ 0 (e^), 


I -VQ 


liQ \-62[\ - vq)^ ^ 

1 /I 7r(2-i/o) 


Po 


+ 0{e^), 

3 - 

8(1-i/o)" ' 2(1-z^o)' 


s + 


+ 0{e^) 


6 

Sijkl = 

n=l 


(5.60) 


for either the spheroid with components of P-tensor (|5.50l] - (j5.52p or the e —>■ 0 limit of the 
spheroid with components of the P-tensor (|5.55p - ()5.59p has components that are related directly 
to the components of the P-tensor via the expressions 

^ I^Piz^o + (1 - z^o)P2^ 

-j ■ 

^ /^2p3Z/o + (1 - Z^o)P4^ 

(-1^27—) • 

S 6 = 2/Z0P6- (5.63) 


Si = 


(pi + 2p2Vo\ 


S 3 — 2/xo 
•S 5 = 2/roP5 


P 3 + Z^ 0 P 4 

1 — 2 z^o 


Finally, it is straightforward, but rather tedious, to show, using the P-tensor derived in the 
previous example, that the strain concentration tensor associated with an isotropic spheroid 
embedded in an isotropic host phase is 

6 

•A-ijM = Cln'K-ijki (5.64) 

n=l 
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where with A = ^qiq^ — (72935 


ai = 


Q!4 = 


94 

2A’ 

91 

2A’ 


and where 


92 

_ 53 

2A’ 

(5.65) 

1 

1 

(5.66) 

O5 = -5 

06 = — 

95 

96 


9i — 1 + “^Pi {^d + Pd) + ‘^P2^d, 
93 = 2p3 (Ad + /Id) + P4Ad, 

95 = 1 + 2p5/ld, 

with 

Ad = Ai — Ao, 


92 = Pi Ad + P2 (Ad + 2/Zd), 

(5.67) 

94 = 1 + 2p3Ad + P4 (Ad + 2/rd) 5 

(5.68) 

96 = 1 -h 2p6/4d, 

(5.69) 

Pd = Pi- Po- 

(5.70) 


The spheroidal cavity result is simply (I5.64ll - (|5.70p with Ai = /ii = 0 and so every occurrence 
of Ad and pd is simply replaced with —Aq and —/xq respectively. The result can be obtained in 
terms of vq alone by using the expression Aq = 2//or'o/(l ~ 2z^o)- 

The average of the concentration tensor over uniform orientations of spheroids can be 
obtained using the result in (IC..5.5h in order to derive an expression of the form 

^ijki T 


5.1.4 Elastic layer 

The result for the spheroid can be used in order to determine the P-tensor and concentration 
tensor for an elastic layer, taking e = 0 in (I5.55p - (I5.59I) . 


Pijkl 


1 f I-2 uq 


2/^0 V 1 “ ^0 


-1/4 I -I/O 

^ijki ' 


6 

ijki I ’ 


Eshelby’s tensor easily follows as 

Sijki — 


^0 qjZ I njA I qj 

^ijkt “T l^idke “T tt 


1 - Z/Q 


^ijk£ 


ijki' 


(5.72) 


(5.73) 


Using the P-tensor the concentration tensor for an isotropic layer is straightforwardly de¬ 
termined as 


•^ijkt 'H-i'ik.f T 


ijki 


Aq — Ai 

Ai -|- 2/xi 


T~iijk£ + 


Ao + 2/io \ .,,4 


Ai -|- 2/ii 


n 


ijki 


4_ "jyS 4_ 

“T ^ijki “T 


(5.74) 


5.1.5 Limiting case of a penny-shaped crack 

For a penny shaped crack, terms up to 0(e) are retained in p5.5511 - ^5.59p to obtain 


Pijki — 


PO 


^(3-4z^o) ^^^i 


TT 

16(1 - vq) 


e 


(Pijki + Pijki) 




/ 1 - 2uo _ 7 r(l -4z^o) A 4 

V2(l-z.o) 8(1-z^o) ; 


+ 


7r(7 - Sz/q) 
32(1 - z/o) 


en 


5 

ijki 


+ 



7r(2-z/o) A 

8(1-1^0) 7 


^ijki 


(5.75) 
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Using this to determine the concentration tensor with /ii = 0 as in the cavity limit, it is 
straightforwardly shown that 


Ajfct — (1 — — -(1 — 


+ 


a/S 

^ijki 


( 41/0(1 - Vq) 1 1^, , n ^ 

7r(l — 21 / 0 ) £ ^ 2 ^ 2f>o)(l - + ^fjH 


+ 


4(1 - i/q) 1 16(3 - i/o)(l - z/q) 


7r(2 - z/q) e 


7r2(2 - z/o)2 




ijk£ 


(5.76) 


where terms of 0(e) have been neglected in (15.761) . This expression has singular behaviour as 
e —)• 0 akin to the potential problem result (I4.4ip and when deriving effective properties for 
distributions of cracks, this singular nature is necessary to yield the correct effective behaviour 
[38| . In fact although 0(1) coefficients have been retained in in (j5.76p only the singular terms 
are required in order to determine effective properties. The expression (I5.76p corrects the 
typographical errors given on p. 104 of |69| . 

A common requirement is the determination of the effective properties of a medium com¬ 
prising penny shaped cracks that are uniformly distributed and uniformly oriented inside some 
host material. Using (IC.55h the associated concentration tensor is shown to be 


37re(l-2z/o) 


8(l-z/o)(5-z/o) 2 

157re(2 - z/q) 


+ 0 ( 1 ). 


(5.77) 


We shall now consider the case of an ellipsoid in an isotropic medium. In order to deal with 
this generally in a tensor setting, ideally an orthotropic tensor basis should be used. Although 
it is possible to write down such a basis, details are rather lengthy and in fact for practical 
computation, it is perhaps most sensible to write down the nine independent components of 
the P-tensor and use matrix computations in the manner described after ^5.1.21 above. 


5.1.6 Ellipsoid in an isotropic host phase 

The nine independent components of the P-tensor for an ellipsoid can be defined in terms of 
the function T(e„;ei,e 2 ) and the semi-axes ratios e„. 

The nine independent components of the Eshelby tensor for an ellipsoid in an isotropic 
medium are usually stated in terms of the four components Sim, 5ii22) •S'nss and 5i2i2 together 
with cyclic properties of the indices, in terms of Imn and Im as defined in (|B.42p - ()B.45p . In 
turn these lead to expressions in terms of the fundamental integral £{en', ei, £ 2 ) via (|B.35p and 
(|B.43P - ()B.46I) . As such use ()5.4p with (I4.2ip and (|B.48I) and employ the properties (IB.43P - 
(jB.46p to derive the following compact forms 


Pun 

Pu22 

A133 

P 1212 


167r/xo(l - z/q) 
1 

167r//o(l - 1 ^ 0 ) 
1 

167r/xo(l - z/q) 
1 

327r/xo(l - z/q) 


/11 + 


1 - 4z/o 


167r/zo(l - z/q) 


h, 


(hi - h), 

(hi-h), 

(I12 + /21) + (Jl + la). 

V 12 W ^ 2 TTfio(l- 120)^ ^ 


(5.78) 

(5.79) 

(5.80) 

(5.81) 


All other non-zero components are obtained by cyclic permutation of the indices in the above 
equations. Those components that cannot be obtained via cyclic permutation are zero, e.g. 
P 1112 = 7^1223 = P 1323 = 0- 

For representations and calculations of the concentration tensor it is convenient to use 
the matrix representation of the tensors. This in discussed in the next §by considering the 
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elliptical cylinder and ribbon crack limits. First however the components of the Eshelby tensor 
are stated, using (II.ip and noting the slightly modified notation for Imn in (|B.42p (i.e. the 
factor of o^) as compared with the standard definition, e.g. Mura [80]. The components are 
expressed as 


Q — ^ 7 1 ^ 7 

871(1-1/0) 871(1-1/0) 

(5.82) 

0 1 1-21^0 , 

871(1 — i/o) ei ^^(1 “ ^0) 

(5.83) 

0 1 1-21^0 , 

871(1 — l/o) £3 871(1 — Z^o) 

(5.84) 

- 16^(1 - + 16^(1 - ^„) + '^'1 

(5.85) 


and permutation rules follow as for the P-tensor. 


5.1.7 Elliptical cylinder and ribbon-crack limit 

In ^4.1.5l it was shown that in the limit as 03 ^ 00 , T(l; 61 , 62 ) —5- 0 and 




a2 

CLl + <12 


1 + e’ 




Qi 

CLl + <l 2 


1 

m’ 


(5.86) 


where e = 02 / 01 . These are used in the expressions for Imn and in Appendix [B] and 
substituted into (|5.78p - (|5.8ip to determine the associated P-tensor components. Since the 
P-tensor is still orthotropic there are nine independent components: 


T’liii = e 


4(l + e)(l-i/o)-(l + 2e) 


P2222 — 


4//o(l - z^o)(l + cf 

4(l + e)(l-z.o)-(2 + e) 


4//o(l - i^o)(l + c-Y 


P1122 — 


P3333 — 0, 
—e 


4/io(l - z^o)(l + cY 
P2233 = 0, 


T 2323 — 


1 


4/io(l + e) ’ 

The Eshelby tensor components follow as 


^1133 = 0 , 

-P1313 = 
-P1212 = 


4 ^ 0(1 + e)’ 

(1 - tio)Y + cY - e 
4/Uo(l -z^o)(l + e)^’ 


'S'liii — e 


/1 -|- 2(1 -|- e)(l — 1 ^ 0 ) \ 

V 2{l-vY{l + eY J’ 


5*3333 — 0, 


53311 = 0, 


52211 


—e -j- 2 i/o(l + e) 

2(l-iyo)(l + eY 


53322 — 0, 


52323 


1 

2(1+ e)’ 


52222 

5i133 

5i122 

52233 

5i313 

5i212 


e -(- 2(1 — i^o)(l + c) 

2(l-i^o)(l + cY 


(l-^o)(l + 6)’ 

—e -|- 2e(l -|- e)r'o 

2(l-iyo)(l + eY 

(l-Ko)(l + e)’ 


2(1+ e)’ 

-e-l- (1 - iio)(l + €Y 

2(1 - iio)(l + eY 


(5.87) 

(5.88) 

(5.89) 

(5.90) 

(5.91) 

(5.92) 


(5.93) 

(5.94) 

(5.95) 

(5.96) 

(5.97) 

(5.98) 


noting that Eshelby’s tensor does not possess the major symmetry, unlike Hill’s tensor. 
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Suppose that the host and elliptical cylinder are both isotropic. In order to determine 
the concentration tensor orthotropic tensors are required. Although it is possible to use an 
orthotropic basis set, it is perhaps most convenient to work with the matrix formulation of the 
tensors and derive the concentration tensor using a symbolic mathematical package such as 
Mathematica. In doing this the matrix formulation [^] of the tensor ^ is employed as noted 
in (I5.42p - (l5.4,'ip . The components of the matrix are so long that to list these here would not 
be beneficial but two very useful limits shall be written down. The elliptical cylindrical cavity 
limit is obtained by taking /xi —)• 0 which yields a matrix form of the tensor (referring to e.g. 
(ICTD l as 


On 

012 

013 

0 

0 

0 

\ 

012 

O22 

023 

0 

0 

0 


013 

023 

O33 

0 

0 

0 


0 

0 

0 

044 

0 

0 


0 

0 

0 

0 

055 

0 


0 

0 

0 

0 

0 

066 

/ 


where 


ail 

ai3 

022 


(1 — t'o)(l + 2e — 2(1 + 6)1^0) 
(1 - 2z^o) 

(2e — 1 + 2(1 — e)iyo)h'o 

(1 - 21 ^ 0 ) ’ 

(1 — i'o)(2 + e — 2(1 + e)z^o) 

e(l - 21 ^ 0 ) ’ 


033 = Ij 


1 + e 
055 — —— 


012 

021 

023 

044 

066 


(1 — uq ){—1 + 2(1 + e)z^o) 
(1 - 2i/o) 

(1 — + 2(1 + e)z^o) 

e(l - 2i/o) 

(2 - e + 2(e - l)z^o)i^o 
e(l - 2z.o) 

1 + e 
2e ’ 

(1 + 6 )^( 1 -^ 0 ) 

2e 


(5.99) 


(5.100) 

(5.101) 

(5.102) 

(5.103) 

(5.104) 


and taking the limit as e ^ 0 yields the ribbon-crack limit, retaining terms up to 0 ( 1 ) in e, 


Oil — 

1 - 

012 = 

013 = 

-z^o, 

021 = 

022 = 

2(1--o)\. .. 

(1 - 2 vo)e ^ 

023 = 

033 = 

1, 

044 = 

055 = 

1 

2’ 

066 = 


-(1 - z^o), 

(5.105) 

2{l-vo)vo . . 

(5.106) 

2(1 — i^o)^o 
(1 - 2z/o)e 

(5.107) 

1 1 
^ 2’ 

(5.108) 

^(1 - I/q) + (1 - l^o). 

(5.109) 


5.1.8 Flat ellipsoid 

Consider the case when oi > 02 S> 03 . It is straightforward to take this limit in (|B.36ll - (IB.40p 
in order to obtain 


11 = 47re2 


{F{k)-E{k) 


h = 4:'K£2E{k) — Ii, 
/a = dvr — 4:'K£2E{k) 


(5.110) 

(5.111) 

(5.112) 


where with reference to (|B.39p and ()B.40p . E{k) and E{k) are introduced as the complete 
Elliptic integrals of the first and second kind, respectively 



dx 

k"^ sin^ x)^0 ’ 



{1 — k^ sin^ x)^/^ dx 


(5.113) 
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and k = 1 — From (I5.110n - (l5.112p . the Imn can be straightforwardly determined via (jB.43p - 
(jB.46p and thus the components of the P-tensor from (I5.78p - (j5.8ip and (j5.82D - (j5.85l] . 

5.1.9 Spheroid limit check 

One can straightforwardly take the spheroidal inhomogeneity limit ai = 02 = a 7 ^ 03 in the 
ellipsoidal result above. In particular it is noted that in the limit as ai —)• 02 = a, referring to 

mM 


where e 


/i = /2 = 27r(l-5(e)). 

(5.114) 

I 3 = 47 r 5 (e), 

(5.115) 

03 / 0 . This then gives 


hi - h 2 - h 2 - L 21 - vr > 

(5.116) 

As - I23- g2 _ 1 ’ 

(5.117) 

% = L - 

(5.118) 

-^31 = h 2 = 

(5.119) 


These can then be used in (I5.78p - (j5.81l] together with the cyclic properties to derive the com¬ 
ponents of the P-tensor for a spheroid. It is straightforward to check that this gives rise to the 
coefficients pi — pe as defined for a TI tensor in (I5.50p - (I5.52D . 


5.2 Anisotropic host phase 

In the potential problem case, scaling coordinate systems assisted in the derivation of results 
associated with anisotropic media. Although such methods can sometimes lead to modest sim¬ 
plifications in elasticity, the general theory does not lead to any significant advances, certainly 
for the problems that are of greatest interest in micromechanics. In particular such methods 
do not lead to significant simplifications for generally transversely isotropic media which is 
a material symmetry of great importance. Therefore to derive the P-tensors associated with 
inhomogeneities in anisotropic host phases, it is best to work with the integral form of the 
P-tensor as defined in (13.91) for an ellipsoid. 

Few explicit results are available in general however since the Green’s tensor cannot gen¬ 
erally be determined analytically. One of the few that can however is that associated with TI 
media. Withers derived the associated Eshelby tensor for an ellipsoid m using the form of 
the Green’s function determined by Pan and Chou [90]. Let us here state his result in the case 
of a spheroid in a TI medium where the semi-major or minor axis of the spheroid is aligned 
with the axis of transverse isotropy. This result shall then be checked by employing the general 
integral form (13.9p . Only in the last decade have articles started to appear that compute effec¬ 
tive properties via micromechanical methods, see e.g. [103] . [32j . It is also important to note 
specihc results for the Eshelby and Hill tensors associated with cracks in anisotropic media. 
See e.g. Gruescu et al. [33| and Barthelemy |3|. 


5.2.1 Spheroid in a transversely isotropic host phase 

Consider a spheroid with semi-axes a = ai = 02 7 ^ 03 embedded in a transversely isotropic 
host phase where the X 1 X 2 plane is the plane of isotropy. The elastic modulus tensor of the 
host is 


^ijki 


E 

n=l 


' ^ijki 


(5.120) 
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where 


c? = 2i^o, cO = 4, c° = 4, (5.121) 

C 4 = no, C 5 = 2mo, Cg = 250 - (5.122) 

Here Kq and mo are the in-plane bulk and shear moduli and go is the antiplane shear modulus 
(often po is used for the anti-plane modulus but this is not employed here in order to avoid 
confusion associated with components pj of the Hill tensor). 


Derivation from Withers’ Eshelby tensor 


Withers derived the Eshelby tensor for a spheroid in a transversely isotropic host medium. In 
order to state this result it is useful to first define the parameters 
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2 ffo)y^' 
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where lo = i'^oiKo + 'mo)Y^‘^- We note that for elastic materials us £ M but vi,V 2 £ C in 
general with V 2 = ui! where an overbar denotes the complex conjugat^ For Uj £ M define 


Vih{vi) = 47rS{vie), 


h{vi) = - 2/3 (ui). 

Vi 


(5.123) 


When Vi £ C, either of the cases in S{vi£) are valid since they are merely an analytic con¬ 
tinuation of the function (of Vi) into the complex Uj-plane. Note that in the case of isotropy, 
4 = Ao = f'o, go = rno = ho, io = no = Ko + mo = ^0 + 2/io and thus vi = V 2 = V 3 = 1. The 
notation /j = /*(!) is therefore appropriate for the isotropic case, as already introduced. 

For a TI host phase dehned by elastic properties (|5.120l) Withers determined the compo¬ 
nents Sijki in the form 


Siiii = ^ [ 250(1 + Mi)vf - mo] LiVihivi) + ^Dmohivs), (5.124) 

i=l 

2 1 

5'ii22 = ^ [ 250(1 + Mi)vi - 3mo] LiVihivi) - -Dmohiv^), (5.125) 

i=l 

2 

5'3333 = 2 ^[4 - noMiv‘f]vlMiLil3{vi), (5.126) 

i=l 

2 

*51133 = 2 ^[4 - noMiv‘f]viLiIi{vi), (5.127) 

i=\ 

2 

*53311 = 2 'Y^[gov‘f{l + Mi) - mo\MiLiV^l3{vi), (5.128) 

i=\ 

1 ^ 1 

*5i313 = ■i^go'^Livfi'^ + Mi){l 3 {vi) - 2MiIi{vi)) + -Dgoh{v3)vl (5.129) 

i=l 


’^This latter point does not appear to have been recognized in the original papers on this subject, e.g. 
\m- An example of a transversely isotropic material for which V 2 = vi € C is zinc with (all in GPa) 
K = 80, £ = 33, n = 50, m = 63, g = 40, for which vi = 1.1284 -|- 0.6465i to 4dp. 
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where 


p 1 {Kq + mo)/vf - go ^ gp - npvf 

^Trgovs^ * 4 + 50 ’ * ^ ^ 87rno5o(^f - ^iK^' 

Note that slightly different notation has been used here from that in |118| and in particular 
the notation has been used whereas |118j used l 2 {vi) for this term in the corresponding 

equations. This is done here to preserve the symmetry with the isotropic case above so that 
as Vi 1 , h{vi) h- 

Via straightforward contraction with the TI compliance tensor Z?o, i.e. the inverse of 
(15.12011 , (jl.ip then yields 


-Pllll + 4122 

4iii — 4 i22 
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?^o(<5'llll + <5 'i122) — 244133 
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2 mo ’ 

■+"o5*3333 — 453311 

A 
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(5.130) 

(5.131) 

(5.132) 

(5.133) 

(5.134) 

(5.135) 


and where A = KqUo — course this calculation could also be done with the help of 

matrices rather than tensor forms. Since the P-tensor is TI however it is rather straightforward 
to write down the TI tensor basis forms 


where 


6 

Pijkl = PnHjikii 

n=l 


6 

Sijki = Y.Snn^M, (5.136) 

n=l 


Pi = 2 P 1111 — P 1212 , P2= Pz = 4i33) (5.137) 

Pi = P 3333 , Pb = 4212, P 6 = 2 P 1313 . (5.138) 


and similarly for the Eshelby tensor with —>■ and Pijki Sijki- 


Derivation from the direct integral form 

As noted above, the P-tensor will itself be transversely isotropic of the form (|5.136l) . Using the 
direct integral formulation of the P-tensor (j3.9ll . let the unit vector ^ pointing to the surface 
of the unit sphere be parametrized by the two angles (p G [0,27r) and i? G [ 0 , 7 r), i.e. 

4 = cosy^sini?, 4 = sfo + sin'i?, 4 = cos'd. (5.139) 


As such, together with (15.120^ . (13.911 becomes 


where 



_ ^ijkl _ 

(1 + (e^ — 1) cos^ d)+2 


sin 'ddipd'd 


^ijke = {Cj^e^ik) 




(5.140) 


(5.141) 


with Nij defined via NikNkj = dij. The components of N are defined as Nij = It 

is straightforward to implement this in a variety of mathematical packages or programming 
languages. 
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Here let us plot the five independent components of the P-tensor, illustrating that the two 
approaches above agree. Let us take the elastic properties to be transversely isotropic and 
choose the material PZT- 7A|. Although realistically this material would normally be chosen 
as the reinforcing phase in a composite, it is appropriate to illustrate the calculations for a real 
material. Its properties are (all stated in GPa) 

K = 121.2, m = 35.8, £ = 73, n = 175, g = 47.2. (5.142) 

In figures [2][3] the components of the P-tensor are plotted, using the explicit form arising 
from the Eshelby tensor and using the direct evaluation of the integral in order to confirm the 
results. 



Figure 2: Plot of the components Pnn (solid black), P 3333 (dashed blue) and P 1313 (dotted 
red) associated with a spheroidal inhomogeneity (varying the aspect ratio s) embedded in 
the transversely isotropic host phase PZT-7A, using the form of the P-tensor derived from 
the explicit form of the Eshelby tensor for this problem. These explicit results are confirmed 
by calculating the components at discrete values of the aspect ratio by evaluating the direct 
integral form of the P-tensor given in ()5.140p . Note that the limiting values as e —>■ 00 
correspond to the circular cylinder result ()5.29p and when e —)• 0 the layer limit is obtained. 



Figure 3: As with Figure [2] but for the components P 1122 (solid black) and F 1133 (dashed red) 
of the P-tensor. 


^This material is Lead Zirconate Titanate, a material frequently used in piezoelectric composites 
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Concentration tensor 

Since the P-tensor is known one can now go on to deduce the associated concentration tensor. 
First assume that the spheroid is isotropic with elastic modulus tensor 

^ijkl = (5.143) 

Since the concentration tensor will be transversely isotropic, it is convenient to write with 
respect to the transversely isotropic tensor basis, i.e. 

6 

Cl,ki = Y.^nKjki (5.144) 

71=1 

where the coefficients are defined in terms of the two independent elastic moduli k,i and 

2 2 2 
c\ = 2 ki + -/XI, c\ = Ki- -/Xi, c\ = Ki- -/XI, (5.145) 

ci = Ki + ^/xi, C5 = 2 /xi, C6 = 2/xi. (5.146) 

Let us employ (|5.12U|1 . together with the form of P-tensor defined in (I5.136p . We shall also 
exploit the properties of the TI basis tensors described in Appendix 1C. 4. 31 (and in par¬ 

ticular the contraction properties in table [2|) , together with the expressions written down in 
(jC.42p - (|C.44h . The inverse of the concentration tensor defined in (j3.23p can then be determined 
in the form 


6 

Ai,ki = Y,^n'H^3ku (5.147) 

n=l 

where upon defining Cn = c\ — 

ai = l+piCi + 2p2Cz, a2 = piC2 + P2CA, da = psci-bp4C3, (5.148) 

d4 = 1 -I- 2p3C2 -FP4C4, d5 = l-hp5C5, d6 = l-bp6 C6. (5.149) 


The tensor Aijki is then inverted, following the procedure in Appendix 1C. 4. 31 to yield the 
concentration tensor 


where 


ai 




2A’ 




6 

Aijki = OinHijkl 

71=1 


(5.150) 


012 = 


02 


05 — 


05 


«3 = (5.151) 

06 = (5.152) 

06 


Alternatively suppose that the inhomogeneity is transversely isotropic with the same sym¬ 
metry axis as the host, i.e. it possesses the elastic modulus tensor of the form (I5.144P but 
where now the constants c)j are dehned in terms of the 5 independent components of this 
tensor. Then the concentration tensor is again defined by p5.150p but of course now with the 
c)j associated with the transversely isotropic cylinder. This indicates the merit of using the 
above notation since one can still use (j5.15nil - (j5.15211 in this case, merely modifying the to 
account for the transverse isotropy of the cylinder. 

As usual, the matrix form of these fourth order tensors can be employed for computational 
efficiency when the problems lack simple symmetries. 
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5.2.2 Circular cylinder in a transversely isotropic host 

Suppose now that the inhomogeneity is a circular cylinder with its cross-section residing in the 
plane of isotropy of the TI host phase. One can arrive at the corresponding P-tensor in two 
ways. The first is to take the limit e = 03/0 —)• oo in the prolate spheroid case in il5.2.1l The 
second way is to recognize that since the anisotropy of the host will not affect the in-plane 
components of the P-tensor, the tensor will simply be the same as that for an isotropic host 
as derived in T5.1.2I but the elastic properties are modified via Aq -|- /tq Kq and /tq —>• rriQ for 
in-plane components and //q 90 for the anti-plane component. Therefore, from (I5.29P 


Pi 


1 

2{Ko + mo) ’ 


P4 = 0 , 


P2 = 0, 

_ Ko + 2mo 
4mo{Ko + mo)’ 


P3 = 0, 



(5.153) 

(5.154) 


The concentration tensor may then be derived by using these coefficients in (I5.148p - (l5.149li 
and the expressions that follow. 


6 Discussion 

6.1 Association with micromechanics 

One of the primary reasons for deriving the Hill or Eshelby tensors and associated concentration 
tensors is to understand a multitude of aspects of the behaviour of inhomogeneous media, 
including their macroscopic constitutive response and so-called effective properties. Following 
a relatively straightforward argument regarding volume averaging, the effective modulus tensor 
C* of an n -|- 1 phase medium with a distinguishable host phase (phase 0) can be stated as 

m, 


n 

C* = C^+ J2MC''-C^)A^ (6.1) 

r=l 

where (pr is the volume fraction of phase r and is the exact concentration tensor associated 
with embedded phase r. This is in contrast to the concentration tensor A. introduced in T3.3I 
which is the concentration tensor associated with an isolated inhomogeneity, i.e. the presence 
of other inhomogeneities is not accounted for in A. As such if the inhomogeneity phases 
are distributed dilutely then one can merely use the approximation A ~ A in (| 6 .ip . Most 
micromechanical methods use a more sophisticated approximation that can account, in an 
approximate manner at least, for interactions. One of the most commonly employed methods 
is the so-called classical self consistent method [H]. Interaction is approximated in this most 
simple self-consistent scheme by taking the host medium in the determination of A to be the 
unknown effective medium. In general then dSI]) gives rise to a nonlinear system of equations for 
the determination of effective properties. In many cases these are not even algebraic equations. 
Furthermore for the self consistent method, one has to make an assumption in advance of the 
symmetry properties of the effective tensor. For example in the case of aligned spheroids C* 
will be TI. 

The textbooks referred to at the end of ^provide an excellent introduction to the numerous 
micromechanical methods, many of which are based on the form of effective modulus tensor 
defined in dSH). A similar form can be deduced for media where there is no distinguishable 
host phase (e.g. polycrystals) and also for media where multiphysics effects are important as 
described in the next section. 

6.2 Beyond the potential problem and elastostatics 

A large number of explicit, compact results associated with the Hill and Eshelby tensors 
for ellipsoidal inhomogeneities, as well as their associated concentration tensors have been 
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collected, stated and in some cases derived. The intention is that this will be of great utility to 
a large number of researchers for implementation in micromechanical and bounding schemes. 
A thorough discussion of both matrix and tensor (where possible due to space limitations) 
formulations has been carried out. Typographical errors in past articles and reviews have been 
corrected and a common notation has been employed. 

Although the general integral forms (13.3p and (13.9p are useful they should generally be 
avoided where explicit forms are available. Gavazzi and Lagoudas [5D] described a numerical 
implementation for elasticity. It should be noted that recently Masson m derived a new form 
of the P-tensor in terms of a single integral, although the integrand is inevitably more complex 
than that in the surface integral in ()3.9p . 

In the literature many of the cases described above are considered as approximations to 
more complicated shaped inhomogeneities. In terms of the derivation of overall effective prop¬ 
erties this is extremely useful, certainly as a first approximation, since it avoids complex com¬ 
putational simulations. However it must be stressed that more advanced analysis is required 
if detailed micromechanical information such as stress concentration calculations close to in¬ 
homogeneities of a complex shape is required [8]. For finite domains, provided the host phase 
is in some sense much larger than the inhomogeneity, if the inhomogeneity is ellipsoidal then 
the temperature gradient field inside the ellipsoid is well approximated as being uniform. The 
inhomogeneity problem associated with bounded domains is described in the book by Li and 
Wang [S5] which summarizes the work in m, 15S]- 

Still remaining in the context of the potential problem and elastostatics, an important 
extension of the inhomogeneity problem is that of the coated inhomogeneity. This problem 
is popular, not least because it arises as a micromechanics problem in the generalized self- 
consistent method (GSCM) [16]. The so-called double inclusion problem dates back many 
decades and was solved approximately by Hori and Nemat-Nasser m although the approx¬ 
imations involved lead to some rather counter-intuitive predictions when used in the GSGM 
m- Exact solutions in the case of concentric spheroids of ellipsoids have been derived by 
Hatta and Taya |36] in the thermal context and Jiang et al. [l2| in two-dimensional elasticity. 
The case of inhomogeneities with radially dependent material properties has been considered 
in Chapter 3 of jUj amongst others. The coated inclusion is also of interest due to its asso¬ 
ciation with the neutral inclusion problem Associated with the coated inhomogeneity is 
the scenario when the interface of an inhomogeneity with the host phase is imperfect [29] , [53] . 
This imperfection can itself be used as the basis for a neutral inclusion m, i- 

It is important to note that when the inhomogeneity becomes very small, i.e. the case of 
a nano-inhomogeneity, then surface energies become non-negligible. This problem has been 
considered by Sharmi and Ganti m and Duan et al. [T8| for example. Including surface 
energies is important in order to incorporate size-dependent effects in effective properties. 
These are absent in classical micromechanical methods that use the standard Eshelby or Hill 
tensors. 

Eshelby’s problem has also been considered in the context of micro-continuum elasticity 
models, which themselves were introduced in order to bridge the gap between continuum and 
atomistic/molecular models |23|. Micropolar (Cosserat) theory has been considered by Gheng 
and He |12j . |13] and Ma and Hu [65] . Micro-stretch theory has been developed by Ma and 
Hu [66]. Strain gradient constitutive behaviour was studied by Gao and Ma [27], [28]. 

The dynamic problem was considered for spheres and cylinders by Mikata and Nemat- 
Nasser [75], m and more generally in m, m- Rate dependence of the Hill and Eshelby 
tensors has been considered by Suvarov and Dvorak [106] and viscoelastic properties have 
been studied by Wang and Weng m by using transform techniques and correspondence 
principles. Nguyen et al. [H] studied cracked viscoelastic solids using the appropriate Eshelby 
tensor. Extensions to plasticity were considered by e.g. [53], [5l], [26] . 

The Newtonian potential and elastostatics problems are canonical problems that can assist 
with the development of coupled (multiphysics) problems. Dunn and Taya [19] . Dunn and 
Wienecke [20] and Mikata [73], m considered the case of piezoelectricity and the prediction of 
the electroelastic moduli. Li and Dunn [56] and Zhang and Soh [121] considered full coupling 
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and the resulting effective moduli associated with piezoelectromagnetic media. The theory 
associated with poroelastic and thermoelastic behaviour was developed by Berryman [5] and 
extended to the anisotropic case by Levin and Alvarez-Tostado [55]. 

Upon closing it should be noted that it is very fortuitous that such elegant and concise 
uniformity results hold for ellipsoidal inhomogeneities. These results allow a large number of 
expressions to be derived analytically and as such the results have been utilized a great deal. 
Having said that there is much work to be done. As has been noted, the Eshelby conjecture is 
still not fully resolved [T], analysis for general shaped inhomogenities continues [3, specifically 
in the context of stress analysis and resulting effective properties and although computational 
methods are powerful, they are still only able to solve elasticity problems for inhomogeneous 
media with an order of 1000 inhomogeneities in “reasonable” times. For use in Monte-Carlo 
schemes this is therefore still computationally expensive. Nonlinear problems in the context 
of finite elasticity still require attention |120j and this applies to coupled problems as well. 
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A Uniform P-tensors for ellipsoidal inhomogeneities 

Fourier transforms can be applied in a straightforward manner to derive forms of the Green’s 
tensors that are useful in the context of deriving properties of the Hill and Eshelby tensors. 
For arbitrary anisotropy in the potential problem the Green’s function takes the form |80j 

G'(z) = [ I - exp(i^| • z) dCdS{^) (A.l) 

levT'* Js2 CijUij 

where the Fourier transform variable with ^ = |^| and where corresponds to ^ = 1, 

i.e. the surface of the unit sphere. Next since 

1 1*00 1*00 1*00 

^ ^ A-2) 

J —oo J —oo J —oo 

the form (jA.ip becomes 

G(z) = ^ / 5(U 2 ) 77 ^^ (^-3) 

O'K Js2 

In elastostatics an entirely analogous approach shows that 

Gni^) = ^[ <5(Uz)iV,,(|)dS(|) (A.4) 

J 

where 

Nik{i)Nkj{i) = dij, = CijkeCjCe. (A.5) 


A.l The potential problem 

Substitute the general form (jA.lli of the free-space Green’s function into (13.211 and it is found 
that the resulting expression mnst be integrated over the inclusion region Vi as well as x E Fi, 
leading to the form 


where 




1 (9^ 

Svr^ dxidxj 



J(^• ^)dS 


(A.6) 


J{p)=[ S{p-^-y) dy. (A.7) 

JyeVi 

Take the simplest case, where Vi is a sphere of radius a. With x E Vi then p < a since ^ is a 
unit vector. It is then recognized that the valne of J(p) is the area of the disc defined by the 
intersection of the plane $ ■ y = p with the sphere Vi. Since \p\ < a, J{p) = 7r(o^ — and 
carrying out the necessary differentiation gives 


Pi 


sphere _ 






4 '^ Js^ 


dS. 


(A.8) 


Since this integral is over the unit sphere 8 “^ and involves only then this shows that the 
P-tensor is independent of x for a spherical inhomogeneity. This was evident from the fact 
that J{p) is a quadratic function of its argnment. In order to perform the integration over the 
unit sphere, introduce the parameters 'd E [0, tt) and E [0, 27r) via 


= cos <7 sin 


^2 = sin (p sin 1 ?, 


^3 = cos?? (A.9) 


and then 
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where 






(A.ll) 


It is straightforward to extend this derivation in order to derive the corresponding result for 
an ellipsoid. Suppose that the ellipsoid is dehned by 1/ = {y : y^(a^a)“^y < 1} where a is the 
second order tensor defined in (|3.5p with aj being the semi-axes of the ellipsoid. With x G lA, 
the only aspect that changes from the calculation for the sphere is that now the function J{p) 
will be the area of the region defined by the intersection of the plane ^-y = p with the ellipsoid 
Vi- It transpires that 


J{P) 


det a 


7r(a^ 



(A.12) 


Importantly the integral is still only over the unit sphere and the result is (|3.3p . 


A.2 Elastostatics 

One can proceed entirely analogously to the transport case in order to derive the following 
representation of the P-tensor for an ellipsoid in the elastostatics context. Using (|A.4D and 
following the same procedure as for the potential problem one obtains ()3.9p . 


B Potential theory 


Two important integrals over ellipsoids arise in potential theory, having important applications 
in micromechanics. Dehne the two functions 




dy, 




- y| dy 


ivi i^-y 

where Vi is the ellipsoidal domain dehned by the inequality 

yl ,yl ^ , 

of o^ 

and of interest is the case when x G Ui. 

Introduce local spherical polar coordinates via yj = Xj + Zj,j = 1, 2, 3 where 


(B.l) 


(B.2) 


zi = rcosifsini}, 


Z 2 = r sinyjsin'd, 


Z 3 = r cos iH, 


(B.3) 


where G [0, vr), G [0, 27r) and r G [0, oo). The surface of the ellipsoid is given by r = p) 

and a differential volume element is dV = sin 'ddrd'ddp so that 


and similarly 
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Dehning the shifted variable ?/; = -d — 7r/2 and upon dehning ri('ip,p) = Ri{'ij; + tt/ 2, p) these 
become 


R7t{2 I*27V 
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At this point note that 

('T^I‘2. r27r 


/ fi^,^)dOd^p= / [f{'tp,(p) + f{ijj,(p + TT)] 

Jo Jo Jo 
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This pairing is useful to argue that certain integrals below are zero. Evaluating the local 
spherical polar coordinates on the surface of the ellipsoid yields 


where 
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noting that A > 0 and C < 0, and thus 
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where the positive root is chosen since — AC > B^ [AC < 0). Hence 


r(x) = 


1 AC - 2B^ + 2HVB2 - AC 


Svr 


'-tt /2 Jo 


A 2 


cos ifdpdif 


(B.9) 

(B.IO) 

(B.ll) 

(B.12) 

(B.13) 

(B.14) 


and 


T(x) = 


IGvr 


.77/2 1-277 g^2^ (J2 

'-77/2 Jo 


A^ 


A3 ^ A 2 


8HVh2 - AC 4BCVB^ - AC 


+ 


^4 I ^3 jcosifdpdif. (B.15) 

The radical contributions to both potentials can be shown to be zero by appealing to (jB. 8 p 
since it transpires that the relevant integrand f{if,p) possesses the symmetry 


f{if,p) = -f{-if,p + 7r), 

The functions T and T thus reduce to 

r(x) = ( 7 T - Ti 

where 
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B.l Closed integral form for r(x) 
Write Fi = Til + ri 2 where 
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X 3 X 1 cos Ip sin V’ cos (/ 7 N cos ip 
OgOf / 


-dpdip. (B.23) 


The contribution from ri 2 is zero - the first term due to 27r periodicity of the integrand in p 
and the second and third terms due to their being odd in ip. Treating T as a function of aj, 
the form of Ti can be exploited, writing 


3 

r(x) = CT - ^ 
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Therefore once T is determined, r(x) straightforwardly follows. As such, introduce A into the 
form of T in (jB.lSp to obtain 
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Next, the evenness of the integrand is exploited in order to write it as 
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Make the substitution sinV' = a^/^/aP^ +t, where t € [0, 00 ) so that 
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where A{t) = 1 /[of + t)(o| + t){a'^ + t). Finally therefore using ()B.24H 
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is a quadratic polynomial in x. It is then found that 
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Setting t = a|s, 'jn = £{^n',£i,^2) which is defined in (I4.22I1 . with = a^/an- Mura (BOj writes 
(jB.29p in the form 


r(x) = ^(-/ + E 44 
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where 
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The link between and £ is then clear: 
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Finally, note that the integrals In (or equivalently ei,£ 2 )) can be expressed in terms 
of elliptic integrals [80] . In particular assuming that oi > 02 > 03 , 
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B.2 Closed integral form for \l'(x) 

One can also derive an integral form for T(x) although such a derivation is rather lengthy. 
Expression (11.38) of Mura |5U] is employed, which establishes that (no sum over i here, with 
sums being shown explicitly for clarity) 


d4'(x) 

dxi 


Svr 


(/-o2/,)-^(/„ 


n=l 



(B.41) 


noting the additional factor of l/( 47 r) here from our modified definition of potentials as com¬ 
pared with Mura. The integral Imn is defined as 


— 27ra^aia2a3 


f 


ds 


s){al +s)A{s) 


(B.42) 


where the slight modification to Mura’s notation should be noted: Imn = cim^mn where I^^ 
is Mura’s definition of this integral as defined in Chapter 11 of [80]. This modification means 
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that expressions are now defined in terms of non-dimensional quantities. In particular it is 
possible to write Imn in terms of Im, In and e^, Cn as follows, 




(,In Im ) 


m ^ n. 


(B.43) 


mn 




Additional relations are noted as 



nm 


(B.44) 


and 



3/i = 3/ii -I- 121 + / 31 , 3/2 = 3/22 + I 12 + Im-, 3/3 = 3/33 -I- /i 3 -I- 123• (B.46) 

Therefore, differentiating ()B.41I1 with respect to Xj,Xk and then xi (in that order) gives 



Note that ()B.47p is however not a fully symmetric fourth order tensor as it should be since 
derivatives should be able to be taken in any order. As such one can enforce in turn major, 
minor then total symmetry m to show that (no sum over repeated coefficients) 




i + j (B.48) 


and odd derivatives are zero. 

C Cartesian coordinates, rotations and tensors 

Cartesian tensors are used throughout this article. Some of their properties are summarized 
shortly, in particular those associated with higher order symmetrized are discussed. Before 
this a brief review of rotations of Cartesian coordinates is given for completeness. 

C.l Rotations of Cartesian coordinate systems 

Consider a hxed Cartesian coordinate system Xi,i = 1,2,3 and an associated Cartesian coor¬ 
dinate system x[,i = 1,2,3 having the same origin, having been rotated arbitrarily in three 
dimensions. The general three dimensional rotation matrix is constructed as a product of 
three rotation matrices, each of which corresponds to a rotation of the axes about a given axis 
in three dimensional space. Begin by rotating anticlockwise about the X 3 axis by use of the 
matrix 



(C.l) 


which generates the rotated coordinate system Xi = Qjj((p)xj. This is followed by an anti¬ 
clockwise rotation about the xi axis by application of the matrix 



1 0 
0 cos 
0 — sin 'i9 I 



(C.2) 
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which generates the rotated coordinate system xi = Q‘lj{'d)Q}j^{ip)xk- Finally an anticlockwise 
rotation about the X 3 axis is performed by application of the matrix 

( cosV’ sin'0 0 \ 

— sinijj cosip 0 I (C-3) 

0 01 / 

so that our required fully rotated system is derived as 

x'i = Qjj{ip)Q‘jk('^)Qiei^)xi 

= Qiii^,'&,ip)x£. (C.4) 

The domains of the Euler angles are 1 } G [0,7r],(^ € [0, 27r),V' G [0, 27r). The angles 1 } and cp 

correspond to the standard angles with identical notation as used in the spherical coordinate 
system defined above. 


C.2 Cartesian tensors in rotated frames 


Employing tensor product notation (g), a second order Cartesian tensor A with components 
Aij can be written 

A = AijOi (g) Oj (C.5) 

where ej is the jth Cartesian unit basis vector. Similarly a fourth order tensor A is defined as 

A = AijkiOi <g) Gj (g) efc (g) e^. (C. 6 ) 


The relation between components of second and fourth order tensors in the rotated system 
x[, generated by application of the general rotation matrix Q to the tensors in the original 
system x,, i.e. A[ - and are 


■^ij — QikQji-'^klj 

-^ijki — QimQ jnQ kpQlq-^mnpqi 


Aij — QkiQijA’j^l, (C-7) 

Aijkl — Qmi QnjQpkQqiAjyinpq. (C-^) 


Often it is useful to determine the average of a general second order tensor, say A over all 
possible rotations of Cartesian axes (uniformly). The natural way to do this is to take a tensor 
with components A'-, diagonal in some coordinate system This frame has been rotated 
from the fixed system Xj. However it is natural to work in a “fixed” coordinate system Xj, 
which can be considered as being obtained from the system x( via a rotation. Indeed we 
have Xi = Qjix'j. The components of AC are diagonal and as such the components Aij are 
dependent on the Euler angles. If one wishes to determine the average of the tensor Aij over 
all such orientations uniformly, it can be done by carrying out the following integration 


87r2 

1 

Stt^ 


rZir pZn rir 

/ / / Aij{ip,'d,ip) sin?? d'dd.Lpdip 

JQ Jo Jo 

r27r r2n /•tt 

/ / / Qkii<f,'&,'ip)Q£jiv’,'J,'J’)A'kf: sin?? dddipdip (C.9) 

Jo Jo Jo 


where the underline denotes orientation averaging. Alternatively there may be some orientation 
distribution function, say p{(p, ??, ip) that weights the importance of certain distributions. A 
weighted orientation average can then be defined as 


I' 27 t i' 27 v rn 

Piy^^^^J^)Qkiip,^,pj)Qejip,'&,'ip)Ake sind dddipdip. (C.IO) 


Note the normalization condition on the weighting distribution 


1 

Svr^ 



p{(p,d,ip) 


sin?? diJdpdip = 1 . 


(C.ll) 


Analogous expressions to (jC.Qp and (IC.lOp hold for the fourth order tensor case of course. 
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C.3 Second order Cartesian tensors 
C.3.1 Isotropy 

The second order identity tensor is lij = 6ij and with a constant, the general second order 
isotropic tensor is therefore Aij = a5ij. Its inverse, with components Aij is Aij = ^5ij. 

C.3.2 Transverse isotropy 

Upon defining the tensor 


©ij — Sij — 5is6j3 (©-12) 

a second order transversely isotropic tensor (with symmetry axis X3) has the form 

Aij — oiOjj + 0i2,5i2,5jii (C.13) 

(C.14) 

The tensor Bjj defined in ()C.12p possesses the following properties: 

©ij — ^ji’i ^ik^kj — ©^j?; ©^t©^J — (C.15) 

Using these properties, the inverse of A has components Aij that can be written 

Aij = —Qij H- Si3^j3- (C.16) 

ot\ 03 

The concept can easily be generalized to other symmetry axes either by use of rotations of 
coordinate axes as in S dC.ip or by use of the notation Uij = rijUj where rij are the components 
of the direction vector associated with the axis of symmetry, so that for the example described 
above = < 5 * 3 . The associated generalization of (IC.12D is therefore 9ij = 6ij — Uij. 

C.3.3 Orthotropy 

A second order orthotropic tensor has the form 

^ij = oii6ii6ji + a2Si2Sj2 + asSisdjs (©-17) 

which has as its inverse 

Aij = — 6ii6ji H- Si2Sj2 H- SisSjs- (C.18) 

CHl OL 2 03 


c.3. 4 Averaging over orientations 

Uniform orientation averaging of second order tensors can be done mechanically via rotation 
tensors as was described above in ^C.ll and written explicitly in (IC.IOI) . Alternatively for simple 
uniform orientation averaging, a simple aspect of tensor analysis associated with invariants can 
be exploited. Averaging uniformly would give rise to an isotropic tensor with components of 
the form 

Aij — 0.5ij. 

where the underline denotes averaging. Performing a contraction in the original tensor gives 
rise to a quantity that does not change with rotations, i.e. a = is an invariant. Therefore 

Akk = a = o5kk = 3a 

so that a = a/3. The components of the averaged tensor therefore take the form 

Aij = ^a6ij. (C.19) 
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If averages need to be taken with respect to some weighting function, then the mechanical 
process of averaging over angles as in (jC.lOD needs to be followed. 

In the case of a transversely isotropic second order tensor Aij with riij = it is easily 

shown that 

Hjj = —ij ~ 3^b'’ (C.20) 

and so taking a uniform orientation average of Aij in ()C.I3D yields 

Aij — 2 (2q:i + as) 5ij. (C.21) 

In the orthotropic case upon taking uniform averages of (IC.lTp it is shown that 

Aij = + ct2 + (C.22) 

C.4 Fourth order Cartesian Tensors 

The following tensors are used extensively in elasticity applications. See Walpole [110] for a 
comprehensive derivation of all associated theory. 


C.4.1 Isotropy 

First define the following tensors 


^ijkt ~ 


(C.23) 

^ijkl 2 T ^ii^jk) '^^ij^kii 


(C.24) 

lijke = 2 i^ik^ji + ^il^jk) = lijki + ^jkt 


(C.25) 

These have the following properties 



rl _ rl r2 _ t2 

^ ij ki ^ kiij f ^ ij k£ ^ kiij f 

lijki — I kiij 

(C.26) 

and 



rl rl _ rl j‘2 r2 — B — f) 

^ijmn^mnkl ■‘■ijki’) ^ijmn^mnkl -‘■ijkii -‘■ijmn-^mnki 

jl _ r| 

^ijmn-^mnki 

(C.27) 


If cr is a second order tensor, whose components are written in the deviatoric/scalar form 

o-jj = cf'ij + ^cr6ij (C.28) 

where a = then 

= (^ij- (C.29) 

Given a fourth order isotropic tensor with components of the form 


Mjki — ^ctilijke + ‘^ct2lijM (C.30) 

then due to (jC.29p 


^ijke(^ke — 3aia6ij + 2a2(T[j. (C.31) 

Introduce a second fourth order isotropic tensor Bijki = ^j^illjkt + then 

AijmnBmnkl ~ Bij^nnAmrikl. ~ (C.32) 

and the inverse of Aijki, denoted by Aijke, such that AijmnAmnke = hjke. is 
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C.4.2 Cubic system 

Define the cubic tensor 


^ijH — + a2lijk£ + Oi^bijkt (C.34) 

and this tensor has the property that 5ijkt = 1 only if i = j = k = i and is zero otherwise. 
Furthermore 


hjmn^mnke — lijkti 


L 


ijmn 


^mnkl ^ijki 


1 

ij ki * 


Writing the inverse of as 


A 


ijk£ — 


al + aius , 1 ,^2 “3 ^ 

2/ I \-^ijk£' ^ijk£ 

a2{cxi-\-as) 0^2 0^2 


(C.35) 


(C.36) 


C.4.3 Transverse isotropy 

We shall use the Hill basis for transversely isotropic (TI) tensors. There are several slight 
variants on this but the Hill basis is used commonly in the micromechanics literature and so 
it appears sensible to adopt it here. This basis set enables a fourth order TI tensor Aijk£ to be 
written in the form 


6 

A,,k£ = Y.an'H%k£ (C.37) 

n=l 

where are constants. We note that in general X 2 7 ^ X^ since contraction of a TI tensor 
with another TI tensor does not result in a tensor with X 2 = X^. The basis tensors are 
defined by 

'^iiki = ®ij^k3^e3, '^ijke = ^kAsbjs, (C.38) 

T^ijke = ^is^js^ks^es, T^ijke = -^^QikQej + QuQkj - Qij&ki), (C.39) 

T^ijke = + Qu^kS^jS + Qjk^£3^i3 + Qj£^k3^i3), (C.40) 

where Qij = 5ij — The notation T-L signifies the Hill basis. 

Let us define the shorthand notation 

n n — rLijpqnpgk£ 

for contraction between the basis tensors defined in ()C.38p - ()C.40p . The contractions defined 
in (IC.4ip are then summarized in table [2l 






nAi 

7^(5) 


~H^ 

~W^ 

~H^ 

0 

0 

0 

0 

7^(2) 

0 

0 


7^(2) 

0 

0 

7^(3) 

7^(3) 

2^(4) 

0 

0 

0 

0 

^(4) 

0 

0 


^(4) 

0 

0 


0 

0 

0 

0 

7^(5) 

0 

■uie) 

0 

0 

0 

0 

0 

7^(6) 


Table 2: The contractions of the basis tensors 
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Note that it is often useful to write the fourth order isotropic identity tensor and basis 
tensors in the Hill TI basis form, i.e. 


^ijki 
^ ij k£ 

r 2 

^ijk£ 


0/6 

^ijk£i 

^ -172 njS I 0-1/4 I 0-1/5 

2 \^ijk£ ^ijk£ ^ijk£ ' ^^ijk£ ^^ijk£ 


T~^ijke + T~^ijke + ^ijke + 



(C.42) 

(C.43) 

(C.44) 


This allows us to define the inverse of the tensor Aijki, Aij^e in a straightforward manner. 
It is 


where 


ai 


04 


04 — 


2A’ 


where A = aiQ! 4/2 — 0 ( 203 . 


6 

Aijki = ^ ^ 

n=l 


02 


O5 


O 2 

1 

as’ 


(C.45) 


03 


06 


as 

~2A’ 

1 

06 


(C.46) 

(C.47) 


C.4.4 Hill’s shorthand notation 

Hill introduced a convenient short-hand notation regarding fourth order tensors. For a fourth 
order isotropic tensor Aijkt. dehned via 

2 

Aijki = (C.48) 

n=l 

Hill denoted it in shorthand notation as H = (oi, 02 ), i.e. the tensor basis is assumed from the 
outset so only the coefficients are to be prescribed. Commonly for linear isotropic elasticity, 
the case when ai = 3 k, 02 = 2/i is assumed where k and /r are as usual the bulk and shear 
moduli. For a fourth order TI tensor Aij^e defined via 

6 

A,,ki = Y.^nH^3H (C-49) 

n=l 

Hill denoted it in shorthand notation as H = (ai, 02 , 0 ^ 3 , 0 ( 4 , 05 , ae). Commonly for linear TI 
elasticity, the case when ai = 2 k, 02 = 013 = i, 04 = 11,03 = 2 m, 03 = 2 p is assumed. 


C.4.5 Orthotropy 

Orthotropic basis tensors are described in the paper by Walpole |107j for example. For practical 
purposes the matrix formulation of fourth order tensors is extremely useful for higher order 
tensors, beginning with orthotropy for example. This formulation is described in 1C.51 below. 
First averaging over orientations is considered. 


C.4.6 Averaging over orientations 

As with the second order case, averages of fourth order tensors can also be taken over all 
rotations. Uniform averaging would give rise to an isotropic tensor with components of the 
form 


—ij k£ 


= alh 


ijk£ 


+ f^^ijke- 


(C.50) 
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Here Ankk = a and Aikik = b are invariants of the original tensor. Therefore A^^kk = a = 3a 
so that a = a/3. Furthermore A^kik = 6 = a + 5/3 so that /3 = 5/5 — a/15. The components of 
the averaged tensor take the form 


A =—4 — (3 A —A 

—ijki ^—ppqq ijki ' \^±l.pqpq ±±ppqqj ^ijki' 

The process of averaging with respect to a weighting function takes the form 


(C.51) 


■y f'‘l7T r27T p7T 

J J J Pi^^^^'^)QrniQnjQpkQqlA!klmn Sm 6 » dddc/dl/ (C.52) 




where arguments on the rotation matrices have been omitted for conciseness. 

A similar approach can be adopted for the uniform orientation averaging of a general fourth 
order TI tensor. Given the general TI form ()C.37p . the averaged tensor will take the isotropic 
form ()C.50h and therefore it must be the case that 


Appqq — 3a — 2q!i + 2 q ;2 + 2 q ;3 + a^, (C.53) 

Apqpq = a 5/3 = CTi + 04 + 2 q ;5 + 2q:6 (C.54) 

which dehnes a and /3 in terms of the coefficients an, i.e. 

^ijkt ~ T 202 T 2 o 3 + a4)l/jk£ T :j^(ai — 2o2 — 203 + 204 + 605 + Qag^lfjki- (C.55) 

o lo 


C.5 Matrix formulation of fourth order tensors 


Matrix representation and manipulation of second order Cartesian tensors is of course trivial. 
It is also often of great utility to represent fourth order tensors in the form of a six by six 
matrices. In particular a general fourth order tensor T with components Tijki with respect to 
a Cartesian basis can be usefully written in matrix form [T] as 


Tllll 

7"i122 

7'i133 

7i123 

7"i131 

7"ii12 

\ 

72211 

72222 

72233 

72223 

72213 

72212 


73311 

73322 

73333 

73323 

73313 

73312 


72311 

72322 

72333 

72323 

72313 

72312 


7i311 

7"i322 

7^1333 

7i323 

7'i313 

7i312 


7i211 

7"i222 

7"i233 

7i223 

7"i213 

7i212 

/ 


(C.56) 


Of particular importance is the ability to carry out the operations of tensor contraction and 
inversion with the matrix form stated here. Defining the matrix [IT] as 


[IT] = 


/ 1 0 0 0 0 0 \ 
0 1 0 0 0 0 
0 0 1 0 0 0 
0 0 0 2 0 0 
0 0 0 0 2 0 
V 000002 / 


(C.57) 


the contraction TijmnTnmke ™ matrix form is [T][IT][T']. Furthermore, defining [T as the 


matrix associated with the inverse of the tensor T it is straightforward to show that [T ^] = 

The matrix [W] can be used in the formulation of matrix forms of the linear elastic consti¬ 
tutive relations, i.e. given the tensor forms 


— CijklOkii 

the equivalent matrix forms as 

b 

C? 

II 

(C.58) 

H = [C][IT][e], 

[e] = [D][W][a] 

(C.59) 
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where [cr] = (fJn a 22 eras a 23 cria ( 712 )^ is the 1x6 column vector of stresses and similarly 
[e] = (eii 622 633 623 ei 3 612 )^ • Multiplying the equations in (|C.59p from the left by [hh] we 
find that 


H = [c][7], b] = [DM 

where 

[C] = [C], m = [W][D][W] = [C]-^ 

and [ 7 ] = [bh][ 6 ] is the engineering strain. 


(C.60) 

(C.61) 
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